EXCESS CHARGE FOR PSEUDO-RELATIVISTS ATOMS IN 
HARTREE-FOCK THEORY 
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Abstract. We prove within the Hartree-Fock theory of pseudo-relativistic 
atoms that the maximal negative ionization charge and the ionization energy 
of an atom remain bounded independently of the nuclear charge Z and the 
fine structure constant a as long as Zct is bounded. 
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A long standing open problem in the mathematical physics literature is the Ion- 
ization conjecture. It can be formulated as follows. Consider atoms with arbitrarily 
large nuclear charge Z, is it true that the radius (see Definition 11.81) and the maxi- 
mal negative ionization remain bounded? A positive answer to this question in the 
non- relativistic Hartree-Fock model has been given by the second author in [2"5] . 
One of the aims of the present paper is to extend the result taking into account 
some relativistic effects. The ionization conjecture for the full Schrodinger theory 
is still open both in the non- relativistic and relativistic case. See [13], [TH], [T7] . 
[6], [7] and |22j for some Z-dependent bounds on the maximal negative ionization. 
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The best result is that N(Z) = Z + 0(Z a ) with a = 47/56 where N{Z) denotes 
the maximal number of electrons a nucleus of charge Z binds (see [5], [7] and |22|L 
As a model for an atom with nuclear charge Z and N electrons we consider (in 
units where h = m = e = 1) the operator 

N 7 1 

H_j>-. (v CArn?W-g!) + v ^ (i) 

i=l 1 11 l<i<j<N 1 J 1 

where a is Sommerfeld's fine structure constant. The operator H acts on a dense 
subset of the N body Hilbert space %f '■= A^ =1 L 2 (M 3 ; C 9 ) of antisymmetric wave 
functions, where q is the number of spin states. The operator H is bounded from 
below on this subspace if Za < 2/tt (see [9] for TV = 1, [5] and [19] for N > 1). 
In this paper we will consider the sub-critical case Za < 2/tt. Let us notice here 
that to define the operator H there is an issue. Indeed for Za < 2/tt the nuclear 
potential is only a small form perturbation of the kinetic energy and hence one 
needs to work with forms to define the operator H. This has been done in detail 
in 0. 

The quantum ground state energy is the infimum of the spectrum of H considered 
as an operator acting on Hi?. In the Hartree-Fock approximation one restricts to 
wave-functions ip which are pure wedge products, also called Slater determinants: 

2p(x 1 ,a 1 ,x 2 ,(J2,...,x N ,a N ) = ^det(ui(x,-,<73))^ =1 , (2) 

with {ui}f =1 orthonormal in L 2 (M 3 ;C«). The Uj's are also called orbitals. Notice 
that ||V'IIl 2 (r 3JV ,c jv ) = 1- The Hartree-Fock ground state energy is 

E HF ( N, Z, a) := inf{q(t/>,^)|^ € Q(H) and ip a Slater determinant}, 

with q the quadratic form defined by H and Q(H) the corresponding form domain. 
One of the main result of the paper is the following. 

Theorem 1.1. Let Z > 1 and a > 0. Let Za — k and assume that < k < 2/tt. 
There is a constant Q > depending only on k such that if N is such that a 
Hartree-Fock minimizer exists then N < Z + Q . 

The idea of the proof is the same as in [53] . One shows that the Thomas- Fermi 
model is a good approximation of the Hartree-Fock model except in the region far 
away from the nucleus. We first introduce some notation in order to introduce the 
Hartree-Fock and Thomas-Fermi models. 

1.1. Notation. Let e be the quadratic form with domain H?(M. 3 ,C q ) such that 
e(u,v) = (E(p)*u,E(p)*v) for all u,v e H^(R 3 ,C q ), (3) 



where E(p) denotes the operator E{iV) — \/ — A + a 2 . As usual (u,v) denotes 
the scalar product of u and v in L 2 (K 3 ,C 9 ). Let V(x) := Za/\x\ and v be the 
quadratic form with domain H 5 (R 3 , C q ) defined by 

v(u,v) = (V?u,vh) for all u, v e H?{R 3 ,C q ). (4) 

From [TO] 5.33 p.307] we have 

I ^^dx<-/ |p||/(p)| 2 dpfor/e^(M 3 ,C) (5) 

JR3 l x l 7T J R 3 

with / the Fourier transform of /. Thus since Za < 2/tt and E(p) > \p\ it follows 
that v(u,u) < e(u,u) for all u e Hi(R 3 ,C q )- 

In the following t denotes the quadratic form associated to the kinetic energy; 
i.e. for all Hi(R 3 ,C q ) 

t(u,v) := a -1 e(u,u) - a~ 2 (u,v) = a" 1 (T(p)^u, T{p)h), (6) 
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with T(p) := E(p) - a~ l . 

A density matrix 7 is a self-adjoint trace class operator that satisfies the operator 
inequality < 7 < Id. A density matrix 7 : L 2 (R 3 ;C 9 ) -> L 2 (R 3 ;C q ) has an 
integral kernel 

7(x,cr,y,r) = ^ XjU^x, a)uj(y, r)*, (7) 

where Aj, Uj are the eigenvalues and corresponding eigenfunctions of 7. We choose 
the Uj's to be orthonormal in L 2 (R 3 ,C 9 ). Let p 7 £ i 1 (R 3 ) denote the 1-particle 
density associated to 7 given by 

/°t( x ) = A jK'( x >°")| 2 - 

cr=l j 

We define 

A := {7 density matrix: Tr[T(p)7] < +00} , (8) 
where for 7 £ .4. written as in Tr[T(p)7] := Tr[E(p)~/} - or 1 Trfy] and 

Tr[£(p) 7 ] :=^A,e( Uj , Mj )- (9) 

Similarly we use the following notation Tr [V7] := ■ Xjv(uj,Uj). 

Remark 1.2. 1/7 G A then p 1 £ L 1 (IR 3 ) since 7 is irace dass and p~ ( e L 4 / 3 (R 3 ). 
Tfte second inclusion follows from Daubechies' inequality, a generalization of the 
Lieb-Thirring inequality (see Theorem \2.3\) . 



1.2. Hartree-Fock theory. In Hartree-Fock theory one considers wave functions 
that are pure wedge products and that satisfy the right statistic: determinantal 
wave functions as in @ . To define the HF-energy functional it is convenient to use 
the one to one correspondence between Slater determinants and projections onto 
finite dimensional subspaces of L 2 (R 3 , C 9 ). Indeed if t/> is given by ([2]) and 7 is the 
projection onto the space spanned by u\, . ..,ujy the energy expectation depends 
only on 7: (ip, Hip) — £ HF (7)- Here £ HF defines the HF-energy functional 

£ HF ( 7 ) = a- 1 Tr[(T(p) - V)j] + X%) - £x (7) , (10) 

where 2? (7) is the direct Coulomb energy 



jk 3 Jr 3 l x _ y| 



and fx (7) is the exchange Coulomb energy 

Jm Jr3 |x - y| 
where we think of the integral kernel 7(x, y) as a q x q matrix. 
Using projections we can define as follows the HF-ground state. 

Definition 1.3 (The HF-ground state). Let Z > be a real number and N > be 
an integer. The HF-ground state energy is 

E HF {N, Z,a) := inf |£ HF ( 7 ) : 7 2 = 7, 7 £ A, Tr[ 7 ] = N) . 

If a minimizer exists we say that the atom has a HF ground state described by 7 HF . 

We may extend the definition of the HF-functional from projections to density 
matrices in A. We first notice that if 7 £ A, then all the terms in £ HF (7) are finite. 
From (|5|) it follows that 



Tr[F 7 ] - ]T A iV ( % , Uj ) < ]T X 3 e( Uj ,u 3 ) = Tr[E(p) 
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On the other hand if 7 G A then p 7 G L^R 3 ) n L% (R 3 ) (see Remark [I~2"]) . By 
Holder's inequality p 7 € L§(R 3 ) and hence ©(7) is bounded by Hardy-Littlewood- 
Sobolev's inequality. The boundness of the exchange term follows from < £ a; (7) < 
2?(7). On the other hand if 7 is a density matrix with j ^ A then £ HF (7) = 00. 
Here we use also that Za < 2/tt. 

Extending the set where we minimize, we could have lowered the ground state 
energy and/or changed the minimizer. That this is not the case follows from Lieb's 
variational principle. 

Theorem 1.4 (Lieb's variational principle, |12j ). For all N non-negative integers 
it holds that 

inf{£ HF ( 7 ) :jeA, 7 2 = 7, Tr[ 7 ] = N} = inf{£ HF ( 7 ) : f € A, Tr[ 7 ] = N}, 

and if the infimum over all density matrices is attained so is the infimum over 
projections. 

The following existence theorem for the HF-minimizcr in the pseudo-relativistic 
case has been recently proved in [2]. 

Theorem 1.5. Let Za < 2/tt and let N > 2 be a positive integer such that N < 
Z+l. 

Then there exists an N -dimensional projection 7 HF = 7 HF (./V, Z, a) minimizing 
the HF- energy functional £ HF given by (JTTJJ) , that is, E HF (N,Z,a) is attained. 
Moreover, one can write 

N 

7 HF ( x >cr,y,T) = ^Ui(x,£7-)«i(y,r)*, 

i=l 

with m € L 2 (R 3 , C), i — 1, . . . , N, orthonormal, such that the HF-orbitals {ui}fL 1 
satisfy: 

(1) h^KFUi — EiUi, with > En > £jv-i > • • • > £1 > —or 1 and 

/i 7 hf := T(p) - ^ + p HF * (xl" 1 - /C 7 hf, (11) 
l x l 

where /9 HF denotes the density of the HF-minimizer and for f G #5(R 3 ) 

N q . 

(£ 7 hf/)(x,<t) = V«,(x,(7) V / Ui(y,r)*/(y,r)|x - y| _1 dy. 

i=l r=l ® 3 

(2) Ui g C°°(R 3 \ {0},C 9 ) /ori = l,...,N; 

(3) Ui G iJ^R 3 \ Bfl(Q)) /or all R>0 and i = 1, . . . , JV. 

In the opposite direction the following result gives an upper bound on the excess 
charge. 

Theorem 1.6. Let aZ < 2. J/iV is a positive integer such that N > 2Z + 1 i/iere 
are no minimizers for the HF-energy functional. 

This theorem for Za < 1/2 was proved by Lieb in jl3j . With an improved 
approximation argument the proof can be extended to Za < 2/n (see [3])- Notice 
that both proofs work not only in the Hartree-Fock approximation but for the 
minimization problem on A Ar L 2 (R 3 ). 

Definition 1.7. Let 7 HF be the HF -minimizer. The function 
^ HF (x):=^- / f!!M dy/or xeR 3 , 

l x l Jn* l x -yl 
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is called the HF-mean field potential and 

$iF(x):= / f!!M dy /or x e M 3 

l x l J\y\<R\K-y\ 

is the HF-screened nuclear potential. 

Definition 1.8. We define the HF-radius F^ F N {v) to the v last electrons by 

f p HF (x)dx = v. 

J\ X \>R*? N (v) 

1.3. A bit of Thomas-Fermi theory. In this subsection we present briefly the 
Thomas- Fermi theory and especially the result that will be used in the rest of the 
paper. We refer the interested reader to [IT] . 

Let U be a potential in L 5 / 2 (R 3 ) + L°°(R 3 ) with 

inf-fllWHoo : U- W € L^(R 3 )} = 0. 
Then the TF-energy functional is defined by 

e?ip) = M s r) % I p^^-l f/(x)p(x)dx + 1 / / ^Mdxdy, 

JR3 J R 3 J r s J R 3 |x - y| 

on non-negative functions p 6 L 5 / 3 (R 3 ) nL 1 (M 3 ). As before, q denotes the number 
of spin states. 

We recall some properties of the TF- model, see [18] , 

Theorem 1.9. Let U be as above. For all N' > there exists a unique non-negative 
P Jj F € L 5 / 3 (R 3 ) such that J pTF < N > and 

£^( P V)=in i {£^(p):peL^(R% f p(x) dx < N>}. 

There exists a unique chemical potential p]j F (N'), with < p]j F (N') < sup U, such 
that p^ F is uniquely characterized by 

= inf{££ F (p)+/^ F (7V') / p(x)dx:0<peL 5 / 3 (R 3 )nL 1 (R 3 )}. 

Moreover p^ F is the unique solution in L 5 / 3 (R 3 ) n L 1 (R 3 ) to the TF-equation 

i(^)l(^ F (x))l = [[/(x) - pl F * lxl" 1 - ^ F (iV')] + . 

If pJj F (N') > then J pJj F — N'. For all fi > there is a unique minimizer 
0<peL 5 / 3 (R 3 )nL 1 (R 3 ) to£^ F (p) + pJ p. 

One defines the TF-mean field potential (pjf , the TF- screened nuclear potential 
^Jj F ji and the TF-radius Rj^ z (i/) to the v last-electron similarly as in Definitions 1 1.71 
and 11.81 replacing the HF-density with the TF-dcnsity. 

Theorem 1.10. IfU(x) — Z/\x\ (the Coulomb potential), then the minimizer of 
£y F , under the condition J p < N, exists for every N. Moreover, p]j F (N) = if 
and only if N > Z . 

When U(x) = Z/\x\ we denote the minimizer of the TF-functional, under the 
condition J p < Z, simply by p TF and J p TF = Z. Correspondingly (p TF and 
<&]j F denote, respectively, its mean field and screened nuclear potential. With this 
notation 

£ TF (p TF ) = -e zI, (12) 
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where eo is the total binding energy of a neutral TF-atom of unit nuclear charge. 

We recall here a result due to Sommerfeld on the asymptotic behavior of the 
TF-mcan field potential, see [23l Th. 4.6]. 

Theorem 1.11 (Sommerfeld asymptotics). Assume that the potential U is contin- 
uous and harmonic for |x| > R and that it satisfies limix^oo t/(x) = 0. 

Consider the corresponding TF-mean field potential tpJj F and assume that /iJj F < 
liminf inf ^ F (x). With ( = (-7 + V73)/2 define 

r~\R |x|— r 

^ F (x) 



a(R) := liminf sup 

r\R | x | =r 



A(R,fi u ) := liminf sup 

r\R | x | =r 



3 4 2- 1 g- 2 7r 2 r- 4 
^ F (x)~^ F 



S^q-^r^ 



- 1 



- 1 



Then we find for all |x| > R 

^ F (x) < ^(l + A( J R, / x7)|x|-C)|x|- 4 + ^ F and 
^f (x) > maxl^Cl + a^lxrO^lxr 4 ,-^)^!- 1 }, 

where 

K^ F ):=inf max{^(l + a(i?)|x|^)- 2 |x|- 3 ,^F| x |}_ 

For easy reference we give here the estimate on the TF-mean field potential 
corresponding to the Coulomb potential. 

Theorem 1.12 (Atomic Sommerfeld estimate, [23, Thm 5.2-5.4]). The atomic 
TF-mean field potential satisfies the bound 

R- min {R'S}^ TF(x) ^ min {^^R}' (13) 

with 2(3 = iri3~i2~^q~i , and for |x| > R > 

^ TF (x)>^(l + a(i?)|xrC r2 | xr 4 , 

where C and a(R) are defined in Theorem ] 1 . 1 11 

Corollary 1.13. Let £ and (3q be defined as in Theorem \ 1 . 1 1\ and \ 1 . 1 2\ resvectivelu. 
Then the TF-mean field potential satisfies the bound 

Z Z% i 
■ if\x\<f3 Zs 

TF / n - J l x | ^PO 



^*(x)> 



3 4 7T 2 



l + aZ-i|x|^)- 2 |x|- 4 if\x.\>p Z- 



2q 2 

with a = /3^(3 2 7r/(g/?|) - 1). 

Corollary 1.14. The TF-screened nuclear potential satisfies 

$TF (x) < 3^d| x |-4 f or aU x g R 3_ 

Corollary 1.15. The following estimate holds 

f (p TF (x))§dx<4||f^zi 
Proof. By the TF-equation and since /x TF = we find 

„TF/„\\| . _ o#/ q \# f /,„TF 



(p"(x))*dx = 2*( g * jr )s / (^-(x))^x. 

JR 3 

The estimate follows from the atomic Sommerfeld upper bound. □ 
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1.4. Construction and main results. We present the basic idea for the proof 
of Theorem 11.11 Let us consider an atomic system with N > 2 fermionic particles 
and a nucleus of charge Z > 1 with Za = k and < k < 2/ir. We assume that 
N > Z and that N is such that a HF-minimizer exists. That is: there exists a 
density matrix 7 HF g .A such that Tr[7 HF ] = N and 

£ HF ( 7 HF ) =inf{£ HF ( 7 ) :7 = 7*,0< 7 </,Tr[ 7 ] = N} . 

Let p TF be the TF -minimizer with potential U(x) — Z/\x\ and under the condition 
J p TF = Z. We know that such a minimizer exist and that the corresponding 
chemical potential is zero (see Theorem II .101) . 

Denoting by p HF the density of the minimizer 7 HF , we find for all r > 

N = f p HF (x)dx 

= f [p HF (x)-p TF (x)]dx+ f p TF (x)dx+/' p HF (x)dx. 

By the equalities above and since /, , <r/ o TF (x)dx < Z, Theorem O follows 
from the following result. 

Theorem 1.16. There exist r > and positive constants c\ and ci independent of 
N and Z but possibly depending on k such that 

f [ ( o HF (x)- ( o TF (x)]dx<ci and f p HF (x)dx < c 2 . 

J\x\<r J\x\>r 

The following theorem is the principal ingredient in the proof of the previous 
one and is the main technical estimate in the paper. 

Theorem 1.17. Let Za = k, < K < 2/vr. Assume N > Z > 1. 

Then there exist universal constants ao > 0, < e < 4 and Cm and C$ depend- 
ing on K such that for all a < ao 

< C $ |x|- 4+E + C M . 

This main estimate is proven by an iterative procedure. We first prove the 
estimate for small x (i.e. |x| < /3o^ _3 ) j then for intermediate x (i.e. up to a fixed 
distance independent of Z) and finally for big x. 

By proving Theorem ll.l7l we also get the following interesting results. The proofs 
of those are given in Section [5] 

Theorem 1.18 (Asymptotic formula for the radius). Let Za = n, < k < 2/ir. 
Both liminf^-foo RWzW) an d um su Pz-s-oo ^WzW) are bounded and behave asymp- 
totically as 

4 22 7T3 _ 1 _1 

3 3 — 2 — v 3 + °{y 3 ) as v ^ co. 
<7 3 

Theorem 1.19 (Bound on the ionization energy of a neutral atom). Let Za = k, 

< k < 2/tt and Z > 1. The ionization energy of a neutral atom i? HF (Z — 1,Z) — 
E HF (Z,Z) is bounded by a universal constant. 

Theorem 1.20 (Potential estimate). Let Za — k, < n < 2/n. For all Z > 1 
and N with N > Z for which a HF minimizer exists with J p HF = N , we have 

|^ TF (x)-^ HF (x)| <A !p \^\-^+A 1) 

with Aq,A\ and Eq universal constants. 



$f x F (x)-$TF (x 



|x| 
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2. Prerequisites 

In this section we recall some results that will be used in the rest of the paper. 
Localization of the kinetic energy. The following is the IMS formula correspond- 
ing to the operator T(p). 

Theorem 2.1 ([H]). Let Xi, i = 0,...,K, be real valued Lipschitz continuous 
functions on K 3 such that X)j=o X? ( x ) = 1 f or a ^ x G R 3 . Then for every f € 

K K 

t(/, /) = tCXi/, Xif) - a- 1 LJ), 
i=0 i=a 
where Li is a bounded operator with kernel 

Mx,y) = ^ Xii ^-_^ r K 2 (a-^-y\), (14) 
where K% is a modified Bessel function of the second kind. 

Remark 2.2. As in [241 App.A, pages 94-98] we use the following integral formula 
for the modified Bessel function 



K 2 (t)=t f 
Jo 



00 

tx/Ff+T 2 



s ds , t > 0. 



We recall that this function is decreasing and smooth in R + . Moreover, 
r+00 

/ t 2 K 2 {t) dt = ?f and K 2 (t) < 16 f" 2 e"5* /or t > 0. (15) 
Jo 

The integral is computed in [211 (A6)] while the estimate follows directly from the 
integral formula for K 2 by estimating \/ s 2 + 1 > \ + | 



,s. 



Generalization of the Lieb-Thirring inequality. This result due to Daubechies 
generalizes the Lieb-Thirring inequality to the pseudo-relativistic case. 

Theorem 2.3 (Daubechies' inequality, [J]). For 7 G .4 

Tr[T(p) 7 ] > f G Q (p 7 (x))dx, 

where G a (p) — ^a~ 4 Cg(a(p/C) 3 ) — a _1 p wi£/i C = .163(7, 9 ^ e number of spin 
states and g(t) = t(l +i 2 )^(l + It 2 ) - ln(< + (1 + t 2 )^). 

Remark 2.4. The function G a defined in the previous theorem is convex and it 
has the following behavior: 



M mui 



in{iaC-ipt,iC-*pt| < Gq ( p ) < § min { I a C7~M , ±CrV } • (16) 



('T/ie proof of the estimate above is in Appendix A.) Notice that when a \, i/ien 
a _1 G Q (p) tends to a constant times p 5 / 3 . 

Theorem 2.5 (Generalization of the Lieb-Thirring inequality, [4]). Lei / _1 oe 

the inverse of the function f(t) := \/t 2 + or 2 — cT 1 , t > 0, and define F(s) = 
J Q dt [/ _1 (i)] 3 . T/ien /or any density matrix 7 it Zioids 



Tr[(T(p) - C/) 7 ] > -Cq / F(|E/(x)|)dx, 

JR3 



witt C < 0.163. 
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Remark 2.6. Since = (t 2 + 2 a - 1 t) 1 / 2 we find for F 

F(s) = 2ia^ 2 ft 3 / 2 (l + \atf' 2 dt for s > 0, (17) 
Jo 

and since by convexity (1 + ^at)i < \[2 + ^(at)i we have 
F{s) < f- a-is§ + ^s 4 for s > 0. 
Hence for any density matrix 7 and potential U G L? (R 3 ) n L A (R 3 ) 

Tr[(T(p) - UH > -CqJ^ (f a~i |C/(x)| * + ^ |C/(x)| 4 )dx. (18) 

Coulomb norm estimate. We present here only the definition of Coulomb norm 
and the result we need. For a more complete presentation we refer to [331 Sec. 9]. 

Definition 2.7. For f,g G Ls(R 3 ) we define the Coulomb inner product 

jr3 it 3 l x - y| 

and the corresponding norm \\g\\c '■= D(g,g)i . 

In the following we write the direct term in the HF-energy functional using 
the Coulomb scalar product: i.e. 2? (7) = Z?(p 7 ,/0 7 ) = D{p 1 ). Similarly, for p G 
L^R 3 ) n ii(R 3 ) the term D(p) denotes D(p,p). 

The next proposition follows as Corollary 9.3 in |23j . 

Proposition 2.8. For s > 0, x G R 3 and f G Li (R 3 ) it holds 

f*M-'< ! [f{y)]J-^---)dy + V2s-^\\f\\ c . 
Moreover, for k > 

/ ^V</ Wkiy + 2h-^-Hf\\c t 

j| y i<ix| i x ~ yi JAQx\,k) i x ~ yi 

where A(\x\,k) denotes the annulus 

A(|x|, fc) := {y G R 3 : (1 - 2fc)|x| < |y| < |x|} . 

2.1. Improved relativistic Lieb-Thirring inequalities. A major difference be- 
tween the pseudo-relativistic HF-model and the non-relativistic one studied in [2 3) 
is that the HF-density p HF in the pseudorelativistic case is not in Ls (R 3 ). By The- 
orem [231 and Remark l2.4l we see that p HF is only in is (R 3 ). Therefore one cannot 
estimate the term p HF * |x| _1 in L 1 -norm simply by Holder's inequality with p = 5/2 
and q = 5/3. To estimate it we are going to use a combined Daubechies-Lieb-Yau 
inequality. 

The following lemma can be found in [24j pages 98-990 
Lemma 2.9. For f G S(R 3 ), 

I £ 1? l/(x)|2dx ^^ L I (/ ' T(p)/) ' 

with p = 7r _1 a -2 . 

The following is a slight generalization of the Daubechies-Lieb-Yau inequality 
formulated in Theorem 2.8 in |24) . 



^The result of the lemma and the proof given in 1241 are actually due to us, but we communi- 
cated the result to the authors of [24] . where it is referred to as a a private communication. 
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Theorem 2.10 (Daubechies-Lieb-Yau inequality). Assume that the potential U € 
Lj oc (R 3 ) satisfies 

> -J7(x) > -kIxI" 1 for |x| < max{a, R} , (19) 
for a, R > and < k < 2/t:. Then we have 

Tr[T(p)-[/]_ > -Cn^a^R^-C^a-i-C f |C/(x)| § + |C/(x) | 4 ) dx. 

J\x\>R 

Proof. If (\/2- 1)/tt < k < 2/vr then K'^a^^R}! 2 + ^a' 1 > Cn^a' 1 and the 
result follows immediately from Theorem 2.8 in [21] observing that for R > a the 
two integrals of the potential on {a < |x| < R} are bounded by the constants. 
If < k < (\/2- I) /it we write 



>R 



U(x) = e-"W [/(^xixKfl + (1 - e" M|a:| )^(i)X|«|<fl + ^ (s)X|x| 
with /i = a~ 2 7r _1 . Using (fT9")l and Lemma [231 we hnd that 

T(p) - C/(x) > iT(p) - «(1 - e-^l^lxl-^ixKfl - tf(x)x H>Jl . 

Hence from the generalization of the Lieb-Thirring inquality Theorem l2.5l fsee (fT5|) ) 
we obtain 



Tr[T(p)-C/]_ > -C [ ^-^(^(l-e-^^l^lxl- 1 )) 1 dx 
-C / ^(l-e-^l^lxr 1 ) 4 ^ 



l\x\<R 

-C [ (a-i|f/(x)|t + |C/(x)| 4 ) dx. 

J|x|>fl 

Since the two first integrals above are estimated below by — Cn 5 ^ 2 a~ 3 ^ 2 R 1 ^ 2 — 
Cn^a" 1 we get the result in the theorem. □ 

By Theorem mU] we find 

k f d y < Tr[T(p) 7 HF ] + C lK Z§Ri + C 2 k 3 Z, (20) 

J\x-y\<R l x _ y| 

with Kg [0, 2/7r], k = Za and i? > parameters to be chosen. This is the inequality 
that we use to estimate p HF * |x| _1 (see proof of Lemma \'S . 2 1 below) . 

2.1.1. Bound on the Hartree-Fock energy. As a first application of Theorem 12.101 
we can give a lower bound to the HF-energy. 

Theorem 2.11 (Bound on the HF-energy). Let N > 0, Z > and such that 
Za = k with < k < 2/ it. Then 

E HF (N : Z) > -2C^Z 2 m - Cn 2 Z 2 ., 
with C the constant in Theorem ] 2.1 (ft 

Proof. Let 7 be a TV-dimensional projection. Since the electron-electron iteraction 
is positive we see that 

£ HF (7) > a -1 Tr[(T(p) - p|)7] 

= a- 1 Tr[(T(p) - ^ X|x|<fl ) 7 ] - a' 1 - X|*|<flh] 

with R > a parameter to be choosen. By Theorem 12. 101 we find 
£ HF (7) > -2C^Z 2 Ni - Ck 2 Z 2 , 
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using that k = Za and by choosing R — C~i Z^ 1 Ni . □ 

3. Near the nucleus 

In this section we prove the estimate in Theorem 11.171 in the region near the 
nucleus (i.e. at distance of Z~s). 

We again assume that N > Z and that an HF-minimizer 7 HF exists for this N 
and Z. We denote the density of 7 HF by p HF . We assume throughout that aZ = n 
is fixed with < n < 2/tt and Z > 1. 

Lemma 3.1. Let Za = k be fixed with < n < 2/tt and Z > 1. Let G a be the 
function defined in Theorem \2.Sl Then, there exists ao > such that for all a < a® 

a' 1 f G a (p HF (x))dx < CZ 7 '\ a- 1 Tr[T(p) 7 HF ] < CZ 7 / 3 

J** (21) 

and \\p TF -p UF f c <CZ 2+ ^, 

with C a universal constant depending only on k. 

Proof. Let p € (0, 1) be such that p~ x n < 2/tt. Notice that here we need k < 2/tt. 
Splitting the kinetic energy into two parts we find 

£HF (7 H F) = ( i_ M)a -iTr[T(p) 7 HF ]+C( 7 HF )-^( 7 HF ) 
+/iTr[(a- 1 T(p)-^) 7 HF ] = ..., 

and introducing p e Li (M 3 ) n L^R 3 ), p > 0, to be chosen 

... = (l-p)a- 1 Tr[T(p)^ F }+p\\p-p iiF \\l ; + (l-p)V( 1 iiF ) (22) 

-£x( 7 HF ) - pD(p) +/iTr[(a- 1 T(p) - (-^ - p* ±)) 7 HF ]. 

p|x| |x| 

Here II • II c denotes the Coulomb norm defined in Definition 12.71 and we used that 



The estimates in the claim will follow from (1221) with different choices of p and 
p. The main idea is to relate, up to lower order term, the last term on the right 
hand side of (|2"2"j) to the TF-energy of a neutral atom of nuclear charge ZpT 1 . This 
has been done in |21| . For completeness and easy reference we repeat the reasoning 
in Propositions IB. II and IB.2I in Appendix [Bj 

To prove the first inequality in (|21l) we choose p as the minimizer of the TF- 
energy functional of a neutral atom with charge pT Z. Since the corresponding 
TF-mean field potential is Z/{p\x\) — p * l/|x| by Proposition IB . 2 1 in Appendix iBl 
we find 

Tr[(a- l T(p) - (-^ - p * ^)) 7 HF ] > -Ctfi + D{p). (23) 

mM l x l 

Here we use $FZ§. Since £ HF ( 7 HF ) < from ^ and fl23J) leaving out the positive 
terms we find 

> (l-M)a- 1 Tr[r(p) 7 HF ]-£a;( 7 HF )-C7 1 Z3. (24) 
From (|2"4")l and Theorem 12. 3l we get 

(l-^a- 1 [ G a (p BF (x))dx<(l-p)a- 1 Tr[T(p)~ f iiF ]<£x(^ F ) + C 1 Z^ (25) 

JR 3 
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It remains to estimate the exchange term. By the exchange inequality (see |15j ) 

£x( 7 HF ) < 1.68 / (^(x))* dx. 

Jr. 3 

To proceed we separate R 3 into two regions. Let us define 

E={xeR 3 :a(C-V HF (x))'>|}, (26) 

with the same notation as in (jTHJ) . By Remark 1 2. 4 [ G a (/? HF (x)) > G2(p HF (x))t in 
£ and a- 1 G Q (p HF (x)) > G 3 (,9 HF (x))§ in R 3 \ E. Hence by Holder's inequality we 
find 

£x( 7 HF ) < 1.68 / (p HF (x))3 dx 

+im({ (^( X ))*dx)*( / p HF (x)dx)^ 

< G 4 / G Q (p HF (x))dx + G 5 ( / a- 1 G Q (p HF (x))dx) I ivi(27) 



Choosing «o such that 1 — fi > 2G4C1 for a < ao, from ((25)) and ((27)) we find 

^a- 1 [ G Q (p HF (x))dx < C x zi +C 5 ( [ a^G^p™ (*)) dxY . 
Jr 3 Wr3 ' 

The first estimate in (|2"Tj) follows from the estimate above using that x 2 — bx — c < 
implies x 2 < b 2 + 2c and that N < 2Z + 1 (Theorem II. 6|) . The second inequality 
in (|2"Tj) follows then from (|23|) and the bound on the exchange term. 

To prove the third inequality in (|21j) we estimate from above and from below 
£ HF (7 HF ). For the one from below we choose in ((22)) p = 1 and p = p TF the 
TF-niinimizer of a neutral atom with nucleus of charge Z. We find 

N 

f HF (7 HF )=E("-(^ lT (P)-^ TF )^) + l^ HF -^ F |lc-^(P TF )-^(7 HF )- (28) 

i=l 

From ((28)) and the proof of Proposition IB. 21 (sec ((6371) 1. we find 

£ HF (7 HF ) > -S.qjddi^^^CZ 2 ^ (29) 
-D(p^) + \\p^-p^\\ 2 c -£ X (^). 

To estimate from above £ HF (7 HF ) we may proceed exactly as in [531 page 543] 
using that a~ 1 T(p) < ^|p| 2 . For completeness we repeat the main ideas. We 
consider 7 the density matrix that acts identically on each of the spin components 
as 

I 3 = / / n p-q d( l d P for 3 = X ' • ■ • > 1- 

33i\p\ 2 <v TF M 

Here n p q is the projection onto the space spanned by /iP ,q (x) :— h s (x ~ q)e 4p x 
where h s is the ground state (normalized in L 2 (R 3 )) for the Dirichlet Laplacian on 
the ball of radius Z~ s with s € (1/3, 2/3) to be chosen. One sees that Tr[7] = Z < 
N since 

p^^^^.hl^^p^.h 2 ^ 

where we have used the TF-equation. Hence £ HF ( 7 ) > £ HF ( 7 HF ). Now we estimate 
from above £ HF (7). Since a~ 1 T(p) < ^|p| 2 and £x{^) > we find 

£ HF (7) < Tr[(-iA - A) 7 ] + D{pj) = . . . , 



EXCESS CHARGE FOR PSEUDO-RELATIVISTIC ATOMS IN HARTREE-FOCK THEORY 13 



and proceeding as in [331 page 543 

- i ff Il«|2 _ zl y2s 



— ii llP \ 2 dpdq-^Z 2s N - I ^p 7 (x)dx + i?(p 7 ). 

Computing the integral and summing and subtracting the term J p TF (/5 TF we get 

£ HF (7) < (^ F (q))*dq-^^ 2s iV- / ^ TF (x)p TF (x)dx 

Jr 3 Jr 3 

-A Pi (k) p TF (x))dx - 2D(p TF ) + D( Pl ). (30) 

R 3 l x l 

By Newton's theorem one sees that D{p 1 ) < D(p TF ) and that 

^W-^W dx < z f p!!M dx < czH^-*). 
m J\x\<z- m 

In the last step we use Holder's inequality and Corollarv ll.151 From (|30|) using the 
TF-equation, that N < 2Z + 1 (Theorem II .6[) and optimizing in s we find 

£ HF (7) < -&,q I (^(d^dd + CZ^-^-D^). (31) 
Jr 3 

Hence from ((29)) and (|3T|) we obtain 

||p HF -p TF ||p<CZ 2 +ft+^( 7 HF ). 

The last estimate in ([2~Tj) follows from the estimate above since £a;(7 HF ) < CZs 
using ([27|) and the estimate just proved on a -1 J G Q (p HF (x)) eix. □ 

Lemma 3.2. Let Za — k be fixed with < k < 2/ir and Z > 1. Then, there exists 
an cto > swc/i i/iai /or a// a < oto, p > and x € R 3 wzt/i |x| < f3Z~~ ^ toe /iaue 



|$IT(x) - $£ F (x)| < a^^(l+ J 53 ^T|x|W+^)| x |- 4+ Tt7i 



Proof. By the definition of screened nuclear potential we have 

*MW *MW ^7 |y|<|x| |x-y| dy -- 
and for all k > by Proposition 12.81 

... < 2i fc -|x|-^H F - P -|| c+ / p HF (y) + ^ F (y) dy , (32) 

° ^A(|x|,fe) |x-y| 
Since ||p TF || L 5 < CZi (Corollary MB and 



/ , 1 |5 rfy <8^|x|^(2fc)i (33) 

JA(\x\,k) |x-y|a 

(see [23] page 549) one finds 

f pTF(y | dy < CZ^|x|5fci (34) 

jA(\ x \,k) l x ~y| 

The term with the HF-density has to be treated differently since we do not have 
a bound for the L 3 -norm of p HF . For a R g R + to be chosen later we consider the 
splitting 

f /° HF (y) f p HF (y) /" p HF (y) 

A( W , fc ) F^ dy = + ^ (l 1<i ^ dy ' (35) 

vl " ' |x-y|>R |x-y|<_R 
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We consider these two terms separately. Let £ be denned as in (j2"6")l ; i.e. the region 
where G Q (p HF ) behaves like (p HF ) 3 (Remark I2.4[) . By Holder's inequality we find 



A(|x|,fe) |x-y| " ~ VM(|x|,*) |x-yl 4 

|x-y|>i?. |x-y|>fl 



yes 



v. I At\~\ t.1 hf — V 2 / V ./../-m>3\ v ' 



U(|x|,fc) |X~y|2 / ^y£R 3 \S 

From the inequality above, Remark 12.41 and estimate ([2~T|) we get 

P HF (y) 



/ 



A(|x|,fe) | X -y 

|x-y|>fl 



■dy < CR-3\x.\*k*Z + C\x\*k*Z?. (36) 



On the other hand for the second term on the right hand side of (|3"5j) by (j2"0)) and 
Lemma 13. II we find 

f f F ( y ] dy < C{zi + R^zi). (37) 

J\x-y\<R l x _ y| 

Hence from (J22), Lemma I3TT1 (|53 |) . <|5B |I and (f3"7 ) . we get 

|$f x F (x)-$T x F ( x )| <C(^j^ + ^|x|*fci +J R-i|x|ifc^ + ^zf +Zf). (38) 

Choosing k such that Z% = Z^\x\^k^, i.e. fc = |x| _1 Z - ^ and i? such that 
.R-sZ 1 -* = zl, i.e. R = Z- 1 we find 

The claim follows using that |x| < fiZ — s . □ 

Theorem 3.3. Let Za = n be fixed with < k < 2/ir and Z > 1. Then there 
exists an tto > such that for all a < ao and x£l 3 with |x| < f3Z~~s we have 

|$£ F (x) - $£ F (x)| < C/3 2 "^(l + f3 2 + /3i + /3 2 +™|x| w)| x |- 4 +m. (39) 
Moreover if |x| < /3.Z /or // < then 

|$£ F (x)-$£ F (x)| < C J 8 a -°('* ) (l + i 9 a +)8*+^W|x| c ^)|x|- ,+0 W, (40) 

OTi/l a(/x) = gg^y - 12 ff M) , &(/*) = 2 + _3_ 24-24m-^ + 1T M flnd C ( M ) = 1_ _ 
.3 — ^49u 

22(8-8^) s ^ r * c "2/ positive constants. 

Proof. Proceeding as in the proof of Lemma 13.21 up to (j36|) we get 

W( x )-*w( x )l ^ C(fc- 1 |xr i ^ 1+ * +2*|x|*ife* +J2-i|x|*feiz) 



|x| W ~|x| 

P HF (y) 
ix- y |<fl l x -y| 



dy, (4i) 



for i? € K + to be chosen. It remains to estimate the last term on the right hand side 
of (|4"Tj) . For 'small' R which is relevant for small x we already did it in Lemma [ 
for 'big' R which is relevant for big x we use Proposition IB. II in Appendix iBl 
Take 7 < 1/263 to be chosen. If |x| < fiZ^^r then by Lemma 1531 



$™(x) - $£ F (x)| < Cp^ (1 + p^^\x\^^)\x\~ 4+ ^ . (42) 
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If instead |x| > fiZ ~ , let H x be the Hamiltonian defined in (|B2|) with P = x 
and v = Z . Then by the definition of H x and taking the HF-minimizer as a trial 
wave function we have 

faf...^,^) < S^i^-Zj ^dy 

ix- y |<H l x - y| 

= M£™( 7 )-z[ fM dy = .... 

^ A y| X - y |<fll x -y| 

Since ^|p| 2 > a _1 T(p), inf 7G ^i £ HF ( 7 ) is estimated from above by the HF-ground 
state energy of the non-relativistic model (i.e. when the kinetic energy is given by 
— 4 A). Moreover, this last one can be estimated from above by £ TF (p TF ) + CW5 Z 2 
(see [IB] and [TJ). Hence we find 

■■■<£ TF (p TF ) + CNiz 2 -Z [ ^\dy. 

J\x-y\<R l x -y| 

On the other hand since |x| > f3Z ~ choosing for some 

from Proposition IB. II it follows that there exists a constant depending only on k 
such that for t G ((1 + 7)/3, min{/, 3/5}), and for every ip € A^L^R 3 ) with 
HV'lb = 1 we have 

(tP,h^) > £™(p™)-c(^ 2 + r 2 )zi-K 

Hence combining the two inequalities above we find 

_HF 



ix-y|<;?. l x - y| 



(43) 



From (|4Tj) and the inequality above we get 

l*w( x )-*w( x )l ^ ck-^i-iz^+czi^ihi 

+CJr*|x|*fc*Z + C{P 1/2 + r 2 )Z^ 3 ' l) . 

Choosing k such that Z^ 3 ~ f ) = zi\x\^ki , i.e k = \x\~ 1 Z^ 1 ~ 5t ^ and R such that 
Zh^-t) ~ R-iz 1+ Te (1 - 5t \ i.e R = fiZ-l+^/A we find 

l*fx|W-^|WI <C(|x|^M* + G8Va + ^-2)^(3-*)). (44) 

Notice that R < f3Z~ l /4 is satisfied choosing / = 4i/3. Then for x such that 
/3Z--HP < | x | < /3Z~i we find 

|cjf x f (x) - $f x F (x)| < C(|x|-S-^S+^ + (^Va + p- 2 )^ 3 ^^^ 3 -^). 

Optimizing in t gives t = 1/3 + 1/99. For this value of t we get 

|$f x F (x) - $£f(x)| < C(l + I3i)l3 2 ~^\ x \-^. (45) 

Inequality (|39p follows from (|4"2"j) and (l4"5"j) choosing 7 such that 4 7 /(l + 7) = 1/66, 
i.e. 7 = 1/263. 

On the other hand from (|44| for x such that j3Z 3 < |x| < j3Z a - we find 
|$f x f(x)-$f x F (x)| < Clxji-^tt+lOjSAlA+l*) 

+c(/? 1 / 2 + /?- 2 )/3^r (^ixr^r (3 - 4) . 

Optimizing in t gives t = 1/3 + 1/99 — j^fx. For this value of t we get 



49m ,1 | 1 49m 



|$f^(x) - (x)| < C(l + pi)(3 2 ~^^ +T ^ JL ^\xr 4+ ^~ - ^7—7 
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Inequality (l4"U)) follows from the one above and (|42p choosing 7 such that 47/(1 

1 49/x 
66(l-£i) 12(1-/j) 



V) - 1 49 M □ 

l) ~ 66(l-u) 12 LJ 



4. The exterior part 

In this section we complete the proof of Theorem 11.171 We first estimate the 
exterior integral of the density and study the minimization problem that the exterior 
part of the minimizer satisfies. Then we prove the main estimate in Theorem ll,17l in 
an intermediate zone, i.e. far from the nucleus but not further than a fixed distance 
independent of Z. To study this area we need first to construct a TF-model that 
gives a good approximation of the HF-density in this intermediate zone. By the 
estimate on the exterior integral of the density we can then also prove Theorem ! 1.1 71 
in the region far away from the nucleus. 

4.1. The exterior integral of the density. The main result of this section is 
the following lemma. 

Lemma 4.1 (The exterior integral of the density). Assume that for some R, a, e' > 


|$f x f(x)-^(x)|<a|xr 4 + £ ', (46) 
holds for |x| < R. Then for < r < R 

(p HF (x) - p TF (x)) dx < ar- 3 +'' (47) 

|x|<r 
and 

f p HF (x)dx < C(l + ar e ')(l + r- s ), (48) 

with C a universal constant. 

We proceed similarly as in the proof of Lemma 10.5 in [53]. Since we need to 
localize we first present some technical lemmas that will take care of the error terms 
due to the localization. The localization error that will appear in the argument 
below (see ([55)1 ) will be in the form of an operator L similar to the error (Ti"4"| in 
the IMS formula. We estimate this error in Lemma [Ol 



Remark 4.2. Let < j3± < .. < ^4 be real numbers with possibly (3^ — 00. Let us 
denote £ r (ft, /3j) = {x e M 3 : for < |x| < fyr}. Then we have 

xes r (ft,ft) K2(a |x y|j rfxdy - 3 fo-fo 
ye £ r (As,&) 

The proof of this estimate is given in AvvendixVA^ 

Lemma 4.3. Let r > and \,v £ (0, 1). Let \- be the characteristic function of 
B r n_i/)(0) and xo be the characteristic function of the sector {x 6 R 3 : r(l — v) < 
|x| < r(l + — A)}. Let r] be a Lipschitz function such thatO < < 1 for all 
x € M 3 , r?(x) = if |x| < r, r?(x) = 1 if |x| > r(l - A) -1 and llVr?^ is bounded. 
Let L denote the operator with integral kernel 

L(x, y) = ^1 (^(^) - ^(y))(^(^)l^l ^ -7(y)lyD ^ 2(c ,-i| x _ yn . m 

4:TT Z |x — y| z 

Then for every function f £ L 2 (R 3 ) we have 

a~ l \(f, Lf)\ < 3D(n, A, r) \\ X of\\l + D(n, A, r)e-\ a ^\\ X -f\\l + cT^/, Qf)\, 
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with D(i],\,r) := \\^v\\oo(j^^zf^~ + an d Q a positive semi-definite operator 
such that 

Tr[Q] < CD(r ] ,X,r)a- 1 r 2 e-h a ' lrv , 
with C depending only on A and v. 

Proof. As a first step we decompose the operator L. We introduce a third cut-off 
function x+ sucn that 1 = x_ (x) + xo ( x ) + X+ ( x ) f° r au x e M 3 . We decompose 
the operator L with respect to these characteristic functions as follows: 

L = X-L(xo + X+) + (xo + X+)Lx- + XqLx+ + X+Lxo + XqLxo- 

We proceed similarly as in [24j Proof of Theorem 2.6 (Localization error)]. For 
Ti,r 2 bounded operators from (ri - T 2 )(ri - T 2 )* > it follows 

TxTl + T 2 Tl < T-lTI + T 2 r|. (50) 

We are going to use several times this inequality with different choices of Ti and 

r 2 . _ _ 

As a first choice we consider Ti = ^fe{x- an d T 2 = l/ v / eT(xo + X+)Lx- with 
E\ > to be chosen. Using ([50)1 we get 

|(/, (x-Hxo + X+) + (Xo + *+)£*-)/) | < ei||X-/lla + -(/, Qif), (^) 

£l 

with Qi — (xo + X+)Lx~L(xo + X+)- We estimate now the trace of Q%. By the 
definition of r),x-iXo an( l X+ it follows that 

T*[Qi]= [ [ L 2 (x,y)dxdy<^^D(r ] ,\,rfr 2 e- a - 1 ™. 

J\x\<r(l-u) J\y\>r 

In the last step we use the definition of L, Remark 14.21 and the definition of the 
constant D(rj, A, r) given in the statement of the lemma. 

Now we choose Ti = y^xo and T 2 = 1/ y/e^X+^Xo with e 2 > to be chosen. 
Proceeding as above we get 

|(/, (x+Lxo + XoL X+ )f)\ < s 2 \\xof\\t + -(/, Qif), (52) 

with Q2 = x+LxqLx+ and such that 

It remains to study the term xo-^Xo- This one has to be treated differently. By 
Schwartz's inequality one gets 

\(f,XoL X of)\<^D(r,,X,r) [ X o(x)|/(x)| 2 , (53) 
since J R3 |L(x,y)|dxdy < ^D(r], A, r). 

_i -1 

The claim follows from ([5T]), and Q choosing e x = D(r],X,r)ae 2" ru \ 
s 2 = ^D(r],\,r) and with Q:=iQi + iQ 2 . □ 

Definition 4.4 (The localization function). Fix < A < 1 and let G : M 3 — s- K be 

given by 

r if |x| < 1, 

G(x):={ f(|x|-l) (1 _ A) 1 _ 1 _ 1 z/l<|x|<(l-A)- 1 , 
I f tf(l-A)- 1 < |x|. 

Let r > and define the outside localization function r (x) := sin(G(^)). 

Remark 4.5. From the definition it follows that ||V0 r ||oo < §"^T^ r_1 - 
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Lemma 4.6. For all r > and A, v € (0, 1) i/ie density p HF o/ £/ie minimizer 
satisfies 

I p HF (x)dx < l + f + 2 sup |x|4> r H ( F _ A) (x)+7^ 

JMyrd-X)- 1 ' |x|=r(l-A) 



'|x|>r(l-A)- 1 |x|=r(l-A) 

with 

K = GD^r- 1 f p HF (x) dx + 2D{\){r- 1 N + Cror 2 )e-\ a ^ rv ', 

with D(X) := (1 + tt/(2A(1 - A)))tt/(2A) and C = C(A, v). 

Proof. Let 7 HF be the minimizer. By the variational principle, 7 HF is a projection 
onto the subspace spanned by u\, . . . ,Un- These functions Ui satisfy the Euler 
Lagrange equations h^vm — SiUi, Si < 0, for i = 1, . . . , N, with h^HF defined in 

Given n a function in C 1 (R 3 ) with support away from zero, we find 

N N 

0>V£ t / Mx)| 2 |x|r? 2 (x)cix = V / u 4 (x)*|x|77 2 (x)/i 7 HFM l (x)dx. 

Since rjT(p)ui 6 L 2 (R 3 ) (Theorem 11.51 (3)), using the Euler-Lagrange equations 
and treating all the terms, except the kinetic energy, as in |23l Formula (63)] we 
get 

N 

> a- 1 VMI • \,rfr(p)ui) - Z / (P F (x)ri 2 (x)dx 

+ / / [p HF (x) P HF (y) - Tr c , | 7 " F (x,y)| 2 ] » = f^fM dxdy 

jr3 ~ i x - y\ 

+ U f p HF (x)7? 2 (x)dx) 2 -± [ p HF (x)^ 2 (x)dx. (54) 

V JR3 7 JR3 

Now we look at the kinetic energy term. For each i £ {1, . . . , N} we may write 

Re(uiV\ ■ \,rfl , (p)u i )=Be(u i v\ ■ I, r(p)(w)) + M^v\ ■ \,[v,T(p)]ui), (55) 

where [A, B] denotes the commutator of the operators A and B. The first term on 
the right hand side of ([55]) is non- negative by the result of Lieb in [T3] . Notice that 
here we may use that r\m £ i? 1 (R 3 ) (see Theorem II .51 (3)). 
Hence, from (l54l) and (l55|) we find 



N 

> a- 1 V Re^ ■ |, [ V , T(p)]m) - Z / p HF (x)r/ 2 (x)dx 

[pH F(x)p H F(y) _ Trc , | 7 H F(x , y) | 2] |y|(l^W(y) ^ 

+ Uj ^(xytxjdx) 2 -!/ p HF (x)^ 2 (x)dx. (56) 

V JS? 7 JR3 

By a density argument we may choose r\ = 8 r the localization function defined 
in Definition 14.41 Reasoning as on page 541 of [53], we get 



AT 

> a 



1 £l^u i » 7 |.|,fo,T(p)]u,) + §( / p HF (x>y 2 (x)dx) 2 

i=l 

+ i+ sup |x|$H ( F _ (x)) / p HF (x)ry 2 (x)dx. (57) 

|x|=r(l-A) 7 JR3 
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It remains to estimate the first term on the right hand side of (|57|) . With the 
same arguments used in the proof of the IMS formula, it can be rewritten as 



A 



iV 



a-^Re^H^TXp)] 



-a 



_1 ^2( u n Lu i), 



(58) 



i=l 



i=l 



where L is the operator defined in (|4"9")l . Using Lemma [4.31 and since ||V?y|| 
|| V^r-Hoo < tt/ (2Ar) we find, with D(X) defined as in the statement, 



N 1 
a- 1 ]T( Ml ,L Ul ) < 3C(A)r- 1 ||xoP HF ||i+^(A)r- 1 e -2^ 1 -||x-p HF ||i 



i=i 



+CD{\)ra- 2 e-^ lr \ 



(59) 



where \q, x- and C are as defined in the statement of Lemma l4.3l Hence combining 
(|57| with (l59|) . using the definition of xo and that Hx-P HF ||i < N we have 



> -3D(A)t 



p HF (x) dx-£)(A)r- x e _ 2 



— 1 — ia x ri/ 



r(l-w)<|x|<r j5 



iV 



-CD(A)ra- 2 e-i Q ' rv + \( ! p HF (x)ry 2 (x)dx) " 

sup |x|d>H ( F _ A) (x)) / p HF (x)ry 2 (x)dx. 

V |x|=r(l-A) ' J«3 

The claim follows using that x 2 — Bx — C < implies x < B + vC- 



□ 



Proof of Lemma We proceed as in [23J page 551]. The first estimate follows 
directly from the equality 

/ (p HF (x) - p TF (x)) dx = ±r f ($* F M ~ <^ F M) du>, 

J\x.\<T JS 2 

and To prove (05]) we use Lemma EH] We first notice that for < (3 < 7 and 
7 such that rj < R 



I P HF (y)rfy < / (p HF (y)-P TF (y))rfy 

J rB<\y\<ri J \y\<r~/ 



I r/3<\y\<r~/ J \y\<ry 

(p HF (y)-P TF (y))dy 

\y\<rff 

< Cr- 3 /3- 3 (l + ar e '). 
Here we used (|4"T)) and that by the TF-equation and (fl"3")) 

„TFf„\ ^ 3 4 2tt 2 



P lt (y)dy 



\y\>r0 



(60) 



|y|>^/3 



Since Ji x i >r P HF _• J| x i>2r/3 '° HF ^° P rove the claim we estimate this second integral. 
By Lemma T4. 61 with r replaced by r/2, A = 4 and ^ = i we get 



|x|>2r/3 



p HF (x)dx < 9 + fr sup $h f /8 (x)+7^, 

|x|=3r/8 



with 7?. defined as in the statement of Lemma 14.61 By (|4"6")l and Corollary 11.141 we 
find 

sup $f /8 (x)<C ff r- 4+£ '+ sup $J F 8 (x)<C(l + ar £ V 4 . 



|x|=3r/8 



|x|=3r/8 
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Moreover, from ([BTJf with (3 — 1/4 and 7 = 1, since N < 2Z + 1 and the boundness 
of R + 3 x H> x p e~ x for all p > 0, we find 

K<C{r- A {l + <rr e ') + r- 1 ), 

The claim follows directly. □ 

4.2. Separating the inside from the outside. We consider the exterior part of 
the minimizer, i.e. the density matrix 

^HF^y^ (61) 

with r as defined in Definition 14.41 This density matrix almost minimizes a new 
energy functional where there is no exchange term. Indeed sufficiently far away 
from the nucleus the electrons are far apart and hence their mutual interaction is 
small. 

We define an auxiliary energy functional on A (see ©) given by 

S A ( 7 ) := Tr[(a-'T(p) - $JP) 7 ] + D( Pl ). (62) 

Theorem 4.7. Let r > and A, v G (0,1). Let xt denote the characteristic 
function of R 3 \ B r (0) . The density matrix 7^? F defined in (1611) satisfies 

£ A {lf F ) < {£ A h) : 7 G Asupp(p 7 ) c M 3 \ B r (0), \\p4t < \\p RF X r\\i} + K, 
where 

K = (^ + § r- l )r~ l f p HF (x) d* + c'a- 2 {l + ar- 2 )e^ a ~ lrd 

+£x( 1 n +C [ [($?(LA)(x)) f + « 3 ($^ A)(X )) 4 1 dx, 

and c' 7 d are positive constants depending only on v and X. 

Proof. We proceed as in J23[ pages 532-6] . The first step of the proof is a localiza- 
tion. Once again we have to treat carefully the localization error coming from the 
kinetic energy. This is the main difference with |23| . For completeness we repeat 
the main ideas of the reasoning. 

We consider the following partition of unity of R 3 : 1 = 9f(x) + 0q(x) + 2 _(x) 
with r defined as in Definition 14.41 and 

0„(x) := (0? (1 _ A) (x) - elWy and 0_(x) : = I 1 - «Ki-A) (*))*• 

Associated to this partition of unity we define 

7o HF ~ #o7 HF #o and 7 HF := ^ 7 HF #_. 

We prove the claim by showing that for all density matrices j £ A such that 
supp(p 7 ) Cl 3 \ B r (0) and ||p 7 ||i < || / o HF x+||i it holds that 

£ A { 1 f F ) + £ HF ( 7 HF ) ~n< £ HF ( 7 HF ) < £ A (l) + £ HF ( 7 HF ). (63) 

The proof of the upper bound in f|63[) is as in [23j page 533] . 

To prove the lower bound as a first step we localize. By Theorem 12. II we find 

N 

a- 1 Tr[T(p) 7 HF ] = a" 1 Tr[T(p)(^ F +^ F )]- a - 1 J>„ (L r + L +£-)«<), 

i=l 

(64) 

where L t ,Lq and L_ are defined as the L,'s in ([13]). 

We first estimate the error term. The procedure is similar to the one used in the 
proof of Lemma |4~51 We introduce three cut-off functions: X- be the characteristic 
function of -B r (i-A)(i-i/) (0), Xr the characteristic function of R 3 \ B r i± JL (0) and Xo 
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defined by Xo( x ) = 1 — Xr(x) — X-( x ) f° r all x & K 3 - Notice that X- and Xr are the 
characteristic functions of sets where 6-, da and 9 r are constants. For k £ {— , 0, r} 
we have the following splitting 

Lk = X-Lk(xo + Xr) + (xo + XrO^fcX- + XrLkXO + XoLkXr + Xo L kXo, 
and proceeding as in the proof of Lemma 14.31 with e\^,£2.k to be chosen we find 

(f,L k f) < e 1 ,k\\x-f\\l+e^ k (f,Qif)+e 2 ,k\\xof\\l + ell(f,Q2f) 

+^l|Vfl*||2o||xo/Hl. 
with operators Q\ and Q2 being positive semi-definite operators with 

Tr[Qi] < g ^)J-») 3 1| VflfcH^e-"-^ 1 -^ 

Tr[Q 2 ] 

Choosing then 



< &^\mr ^ - 



-,r e 



£2,k 



¥1 



and s lik = a \\V6 k kLe 2 



since (||V0 r ||«, + ||V0o||» + l|V0_||») < 37r 2 /(4A>- 2 and ||p HF X-||i < N we get 

a-^^^L. + Lo + i-K) < 3 Sr-V F Xo\\i + 3 Sr- 2 e-^ lrv(1 - X) N 
1=1 

1 -2 -ia _1 ri/(l-A) 

+ ca e 2 v 

Here c is a constant that depends only on v and A. 

Hence from (|6~4"|) , the inequality above and since N < 2Z + 1 we find 

f HF (7 H F) > Tr [( Q -l T(p) _^) (7 HF + 7 HF + 7 HF ) ] +23(7 H F) 



-£x( 7 HF ) - g r- 2 ||p HF xo||i - c'a" 2 (l + «r- 2 )e-2° 
The constants c',d depend only on A and v. Proceeding as in [53] we get 



£ HF ( 7 HF ) > £ HF ( 7 HF )+^( 7r HF )~f 2 :( 7r HF )- C 'a- 2 (l + a r- 2 )e~^" 1 '- d 



Tr 



a- A r( P )-^_ x) (.))^ 



,HF 



^2A x 



37r" „-l\„-l 



2A 1 AX 1 



|x|>r(l-A)(l-f) 



p HF (x) dx. 



The claim follows using Theorem 12.51 



□ 



4.3. Comparing with an Outside Thomas Fermi. At this point we introduce 
an "Outside Thomas Fermi" : a TF-energy functional whose minimizer approxi- 
mates the HF-density at a certain distance from the nucleus. 
Let r > such that 

for all |x| < r for some a > and e' > 0. Let V r be the potential defined by 



K(x)=x+W (x 



if |x| < r, 

.HF/ 



(66) 



<f>™ (x) if |x| > r. 

Here and in the following x^~( x ) := l~Xr(x), x £ R 3 , where Xr is the characteristic 
function of the ball of radius r centered at 0. Notice that V r £ L?(R 3 ) + i°°(E 3 ) 



with 



if{||W||oo :V r -W £ L^(R 3 )} = 0. 
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Let £® TF be the TF-functional £ y F corresponding to the potential V r defined in 
(1661) . Let p® TF be the unique minimizer of £® TF under the condition 



P HF (y)rfy, 

f\>r 



[ p(x)dx< f 

JR3 J\y\ 

(see Theorem ll.9l) . Then pp TF is solution to the OTF-equation 

K^) f (P? TF ) f = [^ TF -M? TF ] + , (67) 



^ TF (x) = K(x)- / ^4dy, 



where 



is the OTF-mean field potential and /4? TF is the corresponding chemical potential. 
From 467} (and p? TF > 0) we see that the support of pp TF is contained in M 3 \ 

B r (0). 

In the intermediary zone instead of comparing directly and $ j^. F we compare 
first the HF-density with the OTF-density and then the OTF-density with the 
TF-density. When comparing the TF and OTF there is no difference with the 
non-relativistic case and for brevity we refer for the proofs to (23j . 

We start by studying the behavior of the minimizer and mean field potential of 
the OTF. The proof of the following bounds is in [231 page 557-558] in the case 
q — 2 and it can be directly generalised to the other values of q. 

Lemma 4.8 ([H Lem.12.1]). For all y £ K 3 we have 

^ TF (y) < 3 4 2-VV|y|- 4 and p TF (y) < S^V^yl" 6 . 



Let (So be as defined in Theorem \l.l'A then for all |y| > /3qZ 3 we have 

^ TF (y)>qy|- 4 andp TF (y)>C\y\- 6 . 
With r,o~,e' such that (|65[) holds and o~r e < 1 we have for all \y\ > r 
P? TF (y) < Cr~ e and ^ r OTF (y) < \V r (y)\ < Cr^. 

Lemma 4.9 ([23l Lem.12.2]). With r,a,e' such that ([65|) holds for all |x| < r we 

have 

(p TF (y)- P BF (y))dy<ar-^'. 

\y\>r 



For x £ K 3 with |x| > r we may write 



(x) - C F (x) = A 1 (r , x) + A 2 (r , x) + A 3 (r , x) , (68) 



where 



and 



A(r,x) = ^-(x)-^-(x), 

A , , f P? TF (y)-p TF (y) , 

A 2 {r,x) = / : : dy 

J\y\>\*\ l x ~yl 



, , , Pp TF (y)-P HF (y) , 

A 3 (r,x.) = / : : dy. 

/r<|y|<|x| l x ^y| 
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4.3.1. Estimate on A\ and A-i. 

Lemma 4.10 ([23, Lem.12.4]). Let N > Z. Given e',a > there exists a constant 
D > such that for all r with (3oZ~3 < r < D for which (f6"5)) holds for all |x| < r, 
then p® TF = and 

^|x|- 4 (l + ar^\^)' 2 < ^ TF (x) < %£|xr 4 (l + Ar<|x|"<) for |x| > r, 
where a, A are universal constants and £ = (—7 + v73)/2. 

Lemma 4.11 ([23, Lem.12.5]). Lei N > Z. Given e',a > iftere ezisis a constant 
D > depending only on e' , a such that for all r with (3qZ~3 < r < D for which 
(|65[) holds for |x| < r, then for all |x| > r 

|Ai(r,x)| < C|x|- 4 - c r c and |A 2 (r,x)| < C\x\- A ~ c r c , 



with £ = (— 7 + V73)/2 and C a universal constant. 

The proof of the previous lemmas is in [23l p. 558-564]. 

4.3.2. Estimate on ||x+P HF - P? TF ||c and Tt[T(p)j^ f } . 

Lemma 4.12. Let G a be the function defined in Theorem \2.3\ and p^ F (x) be the 
one-particle density of the density matrix 7^? F defined in (|61[) . Let Za — n fixed, 
< k < 2/tt and Z > 1. 

Given constants e',a > there exists D < | swc/i t/iai /or aH r wzi/i /3qZ~3 < 
r < D for which (165p holds for |x| < r, it follows that 

a- 1 [ G Q (^ F (x)) dx < a" 1 Tr[T(p) 7r HF ] < 2TZ + Cr~ 7 + Cr~* [ p? F « <*x, 

,/R 3 JR 3 

wii/i C a universal positive constant and TZ as defined in Theorem \4-7\ 



Proof. The first inequality follows directly from Theorem l2.3l To prove the second 
inequality we proceed as in Lemma 13.11 In this case we are interested only in 
the exterior part of the minimizer. Hence, instead of considering the HF-energy 
functional we consider the auxiliary functional £ A , defined in (|62l) . applied to the 
"exterior part of the minimizer" 7 ^ F . 

Splitting the kinetic energy in two terms we find 

£ A (l? F ) > K 1 Tr[T(p) 7j HF ] + D(pf F ) + i TrKa-^p) - 24>« F ) 7r HF ]. (69) 
Since $J? F (x) is harmonic for |x| > r and going to zero at infinity 
$" F (x) < t^j sup $? F (y) for |x| > r. 

|X| |y| =J . 

Hence, since supp(/^ F ) C R 3 \ B r (0) we find 

Tr[(a-'T( P ) - 2d>? F ) 7r HF ] > ^[(^(p) - ^ sup $« F (y))7, HF ] 

I ' I |y|=r 

Adding and subtracting 2D(p,pf F ) for p e L 1 (R 3 ) n ii(R 3 ), p > 0, to be chosen 

... = - v p)l f F ] f [ p " F(x)P i y) <My- (70) 

Jr 3 Jr 3 |x - y| 



where for simplicity of notation here and in the following V p is defined as V p (x) 

l 



^jsup| y | =r $? F (y) - P 



24 ANNA DALL'ACQUA, JAN PHILIP SOLOVEJ 

From ([70]) . (p5)l and the definition of the Coulomb norm and scalar product 
(Definition 12.71) we find 

£ A (lf F ) > \a-^[T{^r] + lD{p™) + ±\\ P ™-p\\l 

-\D{p) + i Tr[(a- l T(p) V P h? F ] (71) 

N 

> Tr[T(p) 7r HF ] + | 2(d r «i, (a- X r(p) - V p )6 rUi ) - \D{p), 

i=i 

denoting by Ui the HF-orbitals. 

We now choose p as the minimizer of the TF-energy functional of a neutral atom 
with Coulomb potential and nuclear charge 2r sup| y | =r <J>^ F (y). Then V p is the 
corresponding TF-mean field potential and we see that the last two terms on the 
right hand side of (I7T1) are like the ones in the claim of Proposition IB. 21 The only 
difference is due to the presence of the localization function 6 r . We now prove that 
these terms give the TF-energy modulo lower order terms. The method is the same 
as that of Proposition [52J We repeat the main steps since in this case the scaling 
depends on r. Notice that since r > /3qZ~3 the contribution is coming only from 
the "outer zone" . 

Let g £ C^°(R 3 ) be spherically symmetric, normalized in L 2 (R 3 ) and with sup- 
port in Bi(0). Let us define g r (x) := r~ 3 g(xr~ 2 ) and tp r := 9r- Since V p is 
subharmonic on |x| > 0, we see from the support properties of Vv and 9 r that 

N N 

^(Mi, (a- x T(p) - V p )6 rUl ) > ^2(e rUi , (a^T(p) - V p * Vv)^) 

i=l »=1 

For p, q £ R 3 we define the coherent states g P ' q (x) := g r (x — q)e* px . By the 



formulas (|B16p and (|B17|) with L q the operator defined in the equation below 
(1517)1 we get 



•■■ = W)*^ 1 I i dpd q (T(p)-aV p ( q ))^2Y,\(0 r ul,g^ 

-a~ l Y, f f dxdq(^)(x)(L q e r «i)(x), (72) 

i=l "' R3 ^ R3 

where u\ denotes the j-th spin component of the orbital it,. By the choice of the 
function g r and with the same arguments that led to (|B19[) in the appendix we find 

a ^ i y2 / dxdq(8 r Ui)(x)(Lq0 r Ui)(x) 

l=1 JR3 JR3 

N 

< 3Y,\\ e r U *\\l\\V9r\\loV0l( S uM9r)) < CV" 4 1 1 P ™ \ | 1 . (73) 



In the first term on the right hand side of (1721 the integrand is zero if |q| < |r 2 
since in this case supp(0 r ) n supp(g q,p ) = (by the choice D < 4/5). To estimate 
it further from below we consider only the negative part of the integrand 



r r N H 

j^a- 1 / dpdq(T(p)-ay p (q))^^|(^,^)| 2 

> ww^ 1 1 1 | q | >ir2 d P rf q ( T (p) - <*v p (<0) , (74) 



T(p)<aV p (q) 
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where we have used that < Yli=i \(^ru{, ffr' q )| 2 — 1 (Bessel's inequality). Wc 
split the domain of integration in p as follows 

{peR 3 : T(p) < aV p (q)} = £1 U E 2 

with Ei,E 2 disjoint and Ei = {p € W. 3 : ^\p\ 2 < V p (q)}. We treat these two 
contributions separately. We have 

p££ 2 P6S 2 

and computing the integral, using that (1 + x)i < 1 + §£C + fa; 2 

■>-C [ dq {a 2 [F p (q)]| + Q 4 [l/ p (q)]}) > -CaV"* -Ca 4 r"¥. (75) 

•/|q|>|r 2 

In the last step we used that [V p (q)] + < 2 j^j su P| x |=r 3>r F ( x ) an d that by the 
hypothesis and Corollary II .141 

r sup $" F (x) < Cr~ 3 , (76) 

|x|=r 

choosing Z? such that ar £ < 1. 

Since T(p) > ±a|p| 2 - ±a 3 |p| 4 we find 



1 Jj M >^ dpdCl ( T (p)-aV p (q)) 



> JJ \*\>\r> dpdq(i|p| 2 -V p (q))-ia 2 yy |q|> , r2 dpdq |p| 4 . (77) 

||p| 2 <V P (q) ^|p| 2 <^ P (q) 

Computing the last integral wc find 

a 2 ff . . 12 dpdq |p| 4 < CarV -1 ^ sup $? F (x))5 < Co?^ . (78) 

Ji, l q l>4 r | x | =r 
^|p| 2 <Vp(q) 

While for the first term on the right hand side of (f77|) , computing the integral with 
respect to p, we get 

lr2 dpdq(i|p| 2 -^(q)) = -4^4 / dq[V P (q)]t 

^|p| 2 <V P (q) lql> ^ 

Hence collecting together §TZ\), ([721), G3) GH> and the inequality above we 

find 

Tr[(a- l T(p)-V p )T?*]>-2l$ [ dx [^(x)]| - Cr- 4 ||/>? F lli " Cr~^ = . . . . 

since (3oZ~^ < r implies fioa^ < n^r. From the TF-equation that p satisfies it 
follows that 

■•■ = U^f) 1 I d*p(*)i-[ p(x)V„(x)dx-Cr- 4 ||^||i-Cr-* 
Jr 3 Jr 3 

= £ TF (p)+D(p) -CV-Wli-CV-*. 

Hence from (fTTj) and the inequality above we get using (IT^I) and (175]) 

f A (7, HF ) > ^-^^(ph^-Cr^-Cr^Hprili. 

The claim follows since £ A {^ F ) < 72- by the result of Theorem 14.71 considering as 
a trial density matrix 7=0. □ 
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Lemma 4.13. Let N' e N and Za = k be fixed, < n < 2/tt and Z > 1. Let 
ej be the first N' negative eigenvalues of the operator a~ 1 T(p) — (p® T F acting on 
functions with support on {x £ R 3 : |x| > r}. 

Given constants e' , a > there exists D < 4/5 such that for all r with (3qZ~z < 
r < D for which (|65p holds for |x| < r, for all [i £ (0, 1) and s < r we have 



]>> > -(Th^T^I [^ OTF (q)]|rfq-Cr-V^-^- 3 r- 5 S 

-(7(1 - aO - ** --8 - (7(1 - n)s- 2 N' : 
with C a positive constant. 

Proof. Let fj be the eigenfunctions (normalized in L 2 (IR 3 ,C 9 )) corresponding to 
the eigenvalues ej, j — 1,..,N'. Let g g Qf(IR 3 ) with support in -Bi(O) and define 
<7s(x) = s~ig(xi/s) for a positive parameter s, s < r. We then write for /i g (0, 1) 

AT' N' 

E e ^ = E(£> («~ lT (p) - ^? TF )/,) = si +^ 

3=1 3=1 

where 

TV' 

si = E^'^-^^p)-^*^)^)' 

3=1 
TV' 

s 2 = Ec/i.^-^-^+^^a/i)- 

3=1 

We estimate these two terms separately. Considering for p,q £ R 3 the coherent 
states .9f' q (x) := e ip x g s (x- q) using (Imp) and (|B17)) . we find 

S i = W //((l-M)«- 1 T(p)-^(q))El(/i.flf' q )| 2 *J* 

3=1 

-(l-^a^W [ dxdqi^xT(L q /,)(x) . (79) 
Estimating the error term as done in (IB32[) and previous inequalities we get 
!: 1 /' ^ E / / rfxdq^x)(L q /,)(x) < (7(1 - m)^ 2 ^V'. 

• =1 JR 3 JR" 



Since we are interested in an estimate from below and (p® TF (q) < for |q| < r, 
from (I79p we find 

B 1 > j^ff ((l-/i)a- 1 r(p)-^ TF (q))El(/3,3 s P ' q )| 2 ^P 

-C(1-m)s _2 ^V'- (80) 
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We estimate now the first term on the right hand side of ([50)1 . Considering only 



the negative part of the integrand and since J2?=i K/j: fff' q )| 2 < 1 we get 



£p [[ ((l-Mja-^pJ-^CqJJ^K/i.^ldq* 
i ( (a ,,^-1tv»^ ,„otf/ 



Now we split the domain of integration in p as follows 

{pe» 3 : a- x (l - /x)T(p) < ^ TF (q)} = E x U E 2) 

with Si,S 2 disjoint and Sj = {p e I 3 : (1 -^)|p| 2 /2 < <^° TF (q)}. We treat these 
two contributions separately. Then 



^ /Xi>n ((1 " ^ a ~ lT ( p ) ~ ^ OTF (q))^q 

pes 2 ' 

// q|> >? TF (q)]+^ = --- 



p6S 2 

> 



(2' 

pes 

and since in the domain of integration 

^[^(0)]+ < IPl 2 < T^[^ F (q)] + (l + ^I)« 2 [^ OTF (q)] + ) 

we get 

£ 2 / > /T,.OTF/ M5 , tfL[,„OTF, 



■■• * "TrVW dq([^Mq)]| + s(i^, (q)]|) 

(WO 2 J|q|>r 

using Lemma \4. 101 in the last step. 

Since V / T+1 2 > 1 + (l/2)i 2 - (l/8)i 4 , we get 

w // q i>,. ((1 - ^) a_lT (p) ~ ^ TF (q))^q 

pes/ 

> /L> r K 1 " ^)ilpl 2 - ^ OTF (q) - K 1 - M)« 2 |p| 4 )rfprfq- 

peSi' 

The last term gives by Lemma [4.101 



|q|» 
pes 



dpdq |p| 4 = a 2 ^ / dq (^) i [^ OTF (q)]| < Ca 2 (^ 

J|q|>r 



'|q|>r 

(82) 

While for the other terms computing the integral with respect to p, we get 

* // q ,>,« 1 -^5lpl 2 -^ F (q))^q = -( T ^) i T5W / ^? TF (q)i 

(83) 

For the term Bi using Theorem 12.51 and Remark 12.61 we find 

B 2 > -CqQT* ||b r OTF - ^? TF * 3 2 ] + llf + aV 3 Hb? TF " ^? TF * 9 2 s ] + \\t). 

2 

From the choice of g s it follows that <^° TF — <^P TF * <? 2 < V r — V r *g 2 and the term 
V r — V r * <? 2 is non-zero only for r — s < |x| < r + s. Hence by Lemma l4.8l and since 
s < r 

||b? TF -^ TF *5 2 ]+lll< / [K(x)-^* 5 2 (x)]^x<Crt, (84) 

2 Jr-s<\x\<r+s 
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and similarly || [(p° TF - </?° TF * 1|| < Cr" 14 s. The claim follows from (|5D|I. (|5T]). 
(155)1 and (gH) using that /3 a^ < n^r. □ 



Lemma 4.14. Let G Q 6e t/ie function defined in Theorem \2.3\ and p^ F (x) the one- 
particle density of the density matrix 7^ F defined in (|6ip . Let Za = k be fixed, 
< K < 2/tt and Z > 1. 

There exists cvq > swc/i i/iai given e',<7 > f/iere exists D < 1/4 swc/i i/iai /or 
aZ/ a < ao <™d r wii/i /3qZ~ 3 < r < D for which (|65l) holds for |x| < r, we have 



\\ x tp™-p™\\ c <Cr-i+* and 
a- 1 



I G Q ( X +p HF (x))dx < Cr-\ a- 1 Tr[T(p) 7r HF ] < Cr~ 7 , 
Js. 3 

with C a universal positive constant. 



(85) 



Proof. The idea of the proof is the same as that of Lemma [XT] In this case we are 
interested only in the exterior part of the minimizer. Hence, instead of considering 
the HF-energy functional we estimate from above and below the auxiliary one £ A , 
defined in (|62p , applied on the "exterior part of the minimizer" 7^ F . 

Step I. Estimate from above on £ A (p/^ F ). Let us consider 7 the density matrix 
that acts identically on the spin components and on each as 



T J = (2¥F / /1 n P>i dpdq ' 

where j € {1, . . . , q} is the spin index, n pq is the projection onto the space spanned 
by ^s' q ( x ) = h s (x — q)e lp x where h s is the ground state for the Dirichlct Laplacian 
on the ball of radius s for < s < r. By the OTF-equation (|S7|) and since p° TF = 
(see Lemma T4.1 Op we see that p 7 (x) = p° TF * |/i s | 2 (x). Moreover, by Lemma \A. 101 

Tr[-IA 7 ] = U^f 3 I (P? TF (x))^ dx + Cs^. (86) 

JR 3 

Since [$? F ]+ G Lf oc (R 3 ), by [H Lemma 8.5] for A' e (0,1) we may find 7 such 
that supp(p 7 ) C {x : |x| > r}, p 7 (x) < /0 7 (x) for xel 3 and 

Tr[(-IA-^ F )7] < Tr[(-iA- X +<f? F )7]+£i / [K(x)]J dx 

+^2Tv) 2 / ArW rfx. (87) 

Since J p 7 < J p 7 = J p^ TF < / X^P HF we may choose 7 as a trial density matrix 
in Theorem 14.71 and we find for A, v to be chosen 

£ A {lf F ) < £ A (l) +n< Tr[(-|A - $« F ) 7 ] + n + D(p^), 

since a~ 1 T(p) < ^|p| 2 . Notice that 1Z depends on A and v. From (|571) it follows 
that 



^(7, HF ) < Tr[(-iA- X +$r)7] + ^i / [K(x)]f dx 

+ Ul^) 2 f Pl (x)dx + n + D( pX( ). 



'W<T 

From the OTF-equation (|67p and Lemma 14.101 we get 



/ p 7 (x) dx < [ P? TF (x) dx < Cr~ 

J\x\<t57 J\x\<^r 
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While since V r (y) < Cr (Lemma 14. 8[) and is non-zero only for |y| > r 



[ [K(x)]| dx < Cr- 



-7 A' 



Hence, from and (|88p and the inequalities above we find choosing A' = r 



£ A {lf F ) < U^P (p? TF (x))tdx- / y r (xK(x)dx + C S ^ 3 

+Cr- 7+ % + K + D(p^) 

Here we used that A' < 1/2 which follows by the bound on D. Since p~ < p 7 , 
D(py) < D{p 1 ). Moreover by Newton's Theorem D(pj) < D(p® TF ). Hence we get 

... < £ OTF (^ TF )+/ K.(x)(p? TF (x)-p 7 (x))dx + C^ 2 r- 3 

+C*r- 7+ i+7e. (89) 

We study now the second term on the right hand side of (|89l) . Since p 1 — p OTF * 
\h s \ 2 , rewriting 

K-(x)(p I OTF (x)- P7 (x))dx= / P ° TF (x)(y r (x)-y r H^I 2 (x))rfx. 



Since s < r, V r is harmonic on |x| > r and p® TF vanishes for |x| < r one sees 
that the integrand on the right hand side of the equation above is non-zero only for 
r < |x| < r + s. Hence by Lemma [ 



K(x)(p^ (x) - p 7 (x)) dx < / p" 1 * (x)K(x) dx < Cr~*s. 

r<\x\<7-+s 



Choosing s = r a we find from (|89|) that 

f A (7, HF ) < f OTF (p? TF ) + Cr- 7 +i + TZ. (90) 

It remains to estimate TL. From Lemma f4. H choosing A, v < 1/2 and D such that 
ar e < 1 we find 



+ p HF (x)dx<Cr- 5 A- 2 . 
•/|x|>r(l-A)(l-*) 

By Lemma 14.81 (|66l) and since A < 1/2 we get 

/ ($H ( F _ A) (x))t dx<Cr- 7 X, 

Jr(l-X)<\x.\< T ^ 



and similarly 



a 3 / ($ r HF _ A) (x)) 4 dx<Cr- 4 A, 
Jr(l-A)<|x|<^ 



since r > (3qZ 3 implies ar 3 < f3 k. Hence from the expression of TZ and the 
boundness of t p e~ t for t > 0, we find 

TZ < £x(lr F ) + Cr- 5 \- 2 +Cr- 7 X. (91) 

We estimate now the exchange term. By the exchange inequality ([TS] or [531 
Th.6.4]) and proceeding as in ([27| we find by Lemma T4. II and Lemma T4. 121 

£z( 7r HF ) < C'[ G Q (p? F (x))dx + Cr-i(a- 1 / G Q (^ F (x))dx) 5 

< CaTZ + Car- 7 + Cr-i(TZ + r- 7 )^. 
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Hence choosing ao such that 1 — Ca > 1/2 for all a < ao we get from the inequality 
above and EU) 



\K < Cr-i{K + r- 7 )^ +Cr- 5 \- 2 + Cr- 7 \, 

that gives 

TZ< C(r- 5 \- 2 + \r- 7 ). (92) 

The second two inequalities in (|55|) follow from the estimate above and lemmas 14.11 
and 14.121 choosing A = 1/2 and replacing r with r/2. 

Step II. Estimate from below on £ A {^ F )- Adding and subtracting D(p® TF ) and 
Tr[p? TF * ^7, HF ] we write 

^( 7j HF ) = TYKa^TCp) - ^ OTF )7, HF ] + IIP? TF - P^Wl - D{p^ F \ (93) 

using that V r — on the support of p^ F . The first term on the right hand side of 
is estimated from below by the sum of the first N' eigenvalues of the operator 
a~ 1 T(p) — (p® TF acting on the functions with support on {x : |x| > r}. Here N' 
denotes the smallest integer bigger than Tr[7^F F ]. Hence by Lemma f4.13l we find for 
p € (0, 1) and s < r 

£ A h? F ) > -(t^) 1 !^/ [^° TP (q)]|rfq-Cr- 8 SM -i-C/x- 3 r- 5 S 

Jm 3 

-C(l - fi)-ir- 5 - C(l - n)s- 2 ( J p? F (x) d X + l) 
+ \\p? TF -P? F \\ 2 c-D(p™) = ... R 

Notice the factor q due to spin. Choosing D such that ar e < 1, by lemmas T4. II and 
l4~T0l we find 

/ p? F (x) dx < Cr~ 3 and / [<^ TF (q)]} dq < Cr- 7 . 
Hence considering p < 1/2 

... > -23 J [cp° TF (q)]ldq-Cr- 7 -Cr- 8 sp-% -Cp- 3 r- 5 s 

-C S - 2 r~ 3 \ ||p r OTF ~ P? F \\ 2 c ~ D(p? TF ) = .... 
By the OTF-equation ([57)1 and since p° TF has support where (p® TF > we find 

■ • • = £ OTF (p r OTF ) - Cr- 7 +i + \\ P r F P? F \\ 2 c, 

choosing p = hr~is~s and s — r 2 ^ . 

Hence combining the inequality above with (|90|) and (|92|) we find 

\\P? TF - P? F Hc < Cr~ 7 ^ + C(r- 5 \- 2 + \r~ 7 ). (94) 
We study now ||x^P HF — P? F ||c- By Hardy-Littlewood-Sobolev inequality we find 

llx+p HF -^ F ||c<qix+p HF -priU<^( / p HF (x)f dx) § . (95) 

To estimate the last term in (|^5]) we are going to use the second estimate in (|55|) 
that we have just proved. With E defined as in (121)1) we find by Holder's inequality 

/ p HF (x)tdx < (f p HF (x)tdx)*(/ Idx 

./ r <| x |< r \yr<|x|, / VJ r <| x |< r 



p HF (x)tdx) M ( / ldJ* 



-<|x|, 

eM 3 \s 
< Cr ~ w A w + C" r - 



r<|x|<^ 
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From the estimate above, and (p?5j) it then follows 

II + HF OTF 1 1 ^ || + HF HP || . || HF OTF n 
\\XrP ~ Pr lie < \\XrP - Pr \\C + \\P r - Pr \\C 

< + C{r- 5 \- 2 + Ar" 7 )2 + C{r~^\^ + r~i\^), 

that gives the claim choosing A = r? □ 
4.3.3. Estimate on A3. 

Lemma 4.15. Let G a be the function defined in Theorem \2.Sl Let Za = n fixed, 
< K < 2/tt and Z > 1. 

There exists ao > such that given e',a > there exists a constant D < 1/4 
depending only on e' and a such that if (|65[) /10/rfs /or a^/ |x| < D, then for all 
a < cto 



a- 1 



f G Q (p HF (y))dy < C|x|- 7 for all |x| < £>, 

J|vl>|x| 



'|y|>l*l 

luii/i C a universal positive constant. 

Proof. If |x| < (3oZ~i we find by Lemma |3~T1 

a- 1 f G a (p HF (y))dy < a" 1 f G a (p HF (y))dy < Czi < C|x|- 7 . 
J|y|>|x| J&3 

While if D > |x| > /3o-Z~^ the claim follows from the second estimate in (|85[). □ 

Lemma 4.16. Le£ Za = k fixed, < ft < 2/-7T, Z > 1 and < /i < ^lg. 

There exists ao sucft i/iai given e' , a > £/iere exists a constant D < 1/4 depend- 
ing only on e' and a such that for all a < ao and for all r with (3qZ s - < r < D 
for which (|65p ZioWs /or |x| < r ; i/ien /or aZZ x with |x| > r 

|^ 3 (r,x)|<c(M)V^, 
zot£/i C > a universal constant. 

Proof. We proceed similarly as in Theorem 13.31 By the formula for ^3 , Proposi- 
tion [575] and Lemma T4. 141 we get 

|-4 3 (r,x)| < / X+Cy) 1 ^^ ~ ^ F(y)l <fr + Cfc-^xl-ir-S+a. (96) 
■/A(|x|,fc) |x-y| 

By Holder's inequality, Lemma T4. 101 the OTF-equation (|67|) and ((33|) we find 

o? TF (y) 



L 



21 . ,1.1 



A(|x|,fc) i x ~ y| 



dy <Cr-^\x\sks. (97) 



Once again, to estimate Jmi x i fc ) Xr ^_ y | ~ ^y we have to proceed differently than 
in [221 Lem.12.7] since /? HF is not in 1/3 (R 3 ). We consider the following splitting 

+< ^p^M, f A), , f p HF (y) , 

Xr(y)i——rdy= ——dy+ ,. __ _d y , 



A(|x|,fe) |x-y J A .(W'*) x-y 7ly|>r, x- y 

|x-y|>fl,|y|>r |x-y|<i? 

(98) 

for i? > to be chosen. By Holder's inequality, Theorem l2.3[ Remark |2~41 (j3"3")l and 
Lemma T4. 141 we get 

' 1 i \ \/ \ 3 3 1121 21 1 1 

dy < Ca^i?"* |x|sfcs r -— + Cr^^\x\^k^ . (99) 



/ A(\x\,k) x-y 

|x-y|>R,|y|>r 
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It remains to study the second term on the right hand side of (|9"5)) . Let v £ R + 
be such that va < 2/tt. We consider the density matrix 7^ defined in (I5T1) with 
A = 1/2. From Theorcm l2.10l it follows that for x such that |x| > r 

TVKcrMp) - ^ XBj?(x) (.)) 7r H / F 2 ] > -C{ylB? + vW). 
Hence we find 

v \ Xr (y I \dy < v \ -dy 

J\y-x\<R l x -y| 7| y -x|<fl|x-y| 

< Tr[a- 1 T(p) 7r H / F] + C(u*R? + v 4 a 2 ) 

and by Lemma T4. 141 

f X^(y) pHF(y ? riy < Cv^r' 7 + C(uiRi + v 3 a 2 ). (100) 

J\ y -*\<R |x-y| 

Hence from <|5B] ) . ([97 ]l . (|99 ]> and (|T00l) it follows that 

|-A 3 (r,x)| < Cv^r' 7 + C(viRi +v 3 a 2 ) + Ca ! iR-i\x\hir-% 



So choosing v = l/2(/3or^ 1 ) (that gives va < 2/tt), k such that r~~ |x| "5 fc = 
/c _1 |x|"2r~2 + 6 : i.e. fc = |x|~T2r 
r _4_ T5|x|T2 i.e. i? = a 2 |x| _ T2r~T 



= 1^ it- i ^ 1- ^ 1 1 5 21 I 1 13 

x i2ri8 and R such that a^R s x s«r 4 + sis 



|^4. 3 (r, x)| < C(r 4+1 -"+|x| 3 ' 6 w-^a + r ^ a 2 + |x|^r 4 ra). 

Finally since r~ 1 a~^ HL < /3 ~ 1 k~s e , the claim follows for |x| > r and fi < 1/(109). 

□ 

4.4. The intermediate region. Here we prove the main estimate in Theorem l 1.171 
up to a fixed distance independent of Z . 

Lemma 4.17 (Iterative step). Let Za = n fixed with < n < 2/tt. Consider 
^ = TT33 an d assume N > Z > 1. 

Then there exists olq > such that for all S,e',a > with S < Sq, where Sq is 
some universal constant, there exists constants £2 , C'^ > depending only on S and 
a constant D = D(e',a) > depending only on e',a with the following property. 
For all a < ao and Rq < D satisfying that j3$Z 3^ < Rq +s and that (|65p holds 
for all |x| < Rq, there exists R' > Rq such that 

\^( X )-^( X )\<C^\-^ 

for all x with Rq < |x| < R' . 

Proof. Let D > depending on a, e' be the smaller of the values of D occurring 
in Lemma 14.111 and Lemma 14.161 Given 6 > 0. We consider Rq < D satisfying 
PaZ- 1 ^ < Rl+ 5 and such that (JBSJ holds for all |x| < R®. 

Set R' = Rl~ s and r = i? 4+<5 . Then we have foZ'i < faZ'^ < r < R < D 
we can therefore apply Lemma 14.111 and Lemma 14.161 From f|68[) we obtain that 
for all |x| > r and all a < ao 

l$f x f (x) - *j2f ( X )| < c|x|- 4 -^ + c( V 4 +^ 

Since for R < |x| < R' we have 

■ 26 r 

<T^< |X|' 



r 
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and thus 

WW - *ff Ml < c|x|- 4 +* + c|xr 4 + 3 T-| xr aoh-*-^). 

Hence choosing Sq sufficiently small there are and £2 such that the claim holds. 

□ 

Lemma 4.18. Let Zat — k fixed with < k < 2/tt. Assume N > Z > 1. 

TTien there exist universal constants ao, £ € (0,4) and D,C$ > 0, D < 1/4, 
such that for all a < ao and x mt/i |x| < D we have 

|$f x f(x)-^(x)|<C $ |x|- 4 + £ . 

Proof. We fix /i = yj as in Lemma 14.171 Since /i < jg , by Theorem 13.31 we 
know that there exists constants a, 6, c > such that for all |x| < f3Z ^ 

|$^(x) - $^( X )| < C(l + /3 2 + /3 5 / 2 + /3 fc |x| c )^ a |xr 4+Q . (101) 

We first show that we may choose 5 small enough such that if we choose R 1+& — 
/3qZ 3- we have for all |x| < R that 

|$f x f(x)-$f x f(x)|<^|x|- 4 +§. (102) 

Let (3 > be such that {[3Z-^r) 1+s = A)^"^, i.e. /3 1+l5 = P Q Z sh P . Hence 
from ([TUT]) we find for all |x| < f3Z-^r 

|$f x J(x) - $£f (x)| < C(l + /? 2 + /? 5/2 + /3 b |x| C )/3 2 " f ^ §i ^|xr 4+ ^ , 
and by the choice of /3 (and (3q < 1) 

|*^(x)-*j^"(x)| < ca + ^^+z^^+z^^+^z-^) 

Hence if <5 is small enough we may choose a universal constant C$ such that (1102[) 
holds. 

Let now 8 be small enough so that we may apply Lemma 14.171 This give constant 
£2 and (depending only on S) and for all <j,e' > a constant D < 1/4. Now 
choose a — max{ C^,C$} and e' — min{a/2,e 2 }. Now a, e' and -D are universal 
constants. To prove the claim we shall prove that for all |x| < D 

|$f x f(x)-$^(x)|<a|xr 4 ^'. (103) 

We have to prove that D belongs to the set 

M = {0 < R < 1/4 : Inequality (ITU5|) holds for all |x| < R}. 

We reason by contradiction. If this was not true then D > Rq = sup.M and in 
particular Ro < 1/4. From (|102[) and the choice of a and e' it follows that either 
i?>l/4ori?eAL In the first case then Rq = sup M. = 1/4 > D that contradicts 
our hypothesis. On the other hand if R £ M, then R^ +s > R 1+s = faZ' 1 ^ . 
It then follows from Lemma [4.171 that there exists R' 6 M. with R' Q > i?o- This 
contradicts also our hypothesis. □ 

4.5. The outer zone and proof of Theorem 11.171 The proof of Theorem 1 1.1 71 
follows directly from Lemma [4.181 and the following result. 

Lemma 4.19. Let Za = k, < k < 2/ir. Assume N > Z > 1. Let D,e and C$ 
be the constants introduced in Lemma \4-18\ 

Then there exist ao > and a universal constant Cm > such that for all 
a < ao and x with |x| > D we have 

|$f x f(x)-$f x f(x)|< C M . 
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Proof. Here Ci, i = 1, . . . , 6 denote positive universal constants. We write 

|$f x f(x)-$TF (x) | < |$HF (x) _ $ TF (x)|+ I P TF ( y)+ Pj^ dy . (104) 

Since $^ F (x) — <j>J, F (x) is harmonic for |x| > D and tends to zero at infinity we 
have by Lemma 0331 

|4>g F (x) - $£ F (x)| < sup |$g F (x) - $£ F (x)| < C^D-^. (105) 

|x|=£> 

For the second term on the right hand side of (| 104[) we write 

P TF (y) + P HF (y) 



£><| y |<|x| l x -y| 



■dy 



|y|>U 

By Lemma 14.11 Lemma 14.181 estimate (fT3"]) and the TF-equation we find 

/ (p TF (y) + p iiF (y))dy<C 1 (l + C^)(l + D- 3 ) + C 1 D- 3 . (107) 

JD<\y\ 

It remains to estimate the first term on the right hand side of (|106D . By Holder's 
inequality, estimate (1131) and the TF-equation we get 

S^ gl &*' s c <Ly' w>l *o ,Di s c = D " 4 - <io8) 

To estimate the term with the HF-dcnsity we use Theorem 12.101 Let 7§ F be the 
exterior HF-density matrix as defined in (|61[) with r = D/2 and A = 1/2. Then by 
Theorem [Uni with v = PlD^ 

a- 1 Tr[(T(p) - -^- X b d (x)(-)hg F 2 ] > -C 4 (£>M + ^ 4 a 2 ), 

|X — •] T 

and thus 

/ ^4 ^ ^ ^a" 1 Tr[T(p) 7 g F 2 ] + C 6 D-\ 

J\*-y\<D/4 |x-y| 7 

Here we use that D > 2f3oZ~i (for a < a>o) and D < 1/4. By Lemma [4.141 we 
conclude 

I xtly) pyl dy <[ f^ dy < c ,n-'. (109) 

J|x-y|<_D/4 l X_ y| J\x-y\<D/4 l x ~ Yl 

The claim follows collecting together formula (1104j) to formula (|109[) . □ 

5. Proofs of Theorems [TTTl fTTTHl H~T$1 and POUl 

In this section we always assume the following: Za = k with < k < 2/tt and 
N > Z > 1. 

Proof of Theorem \l.l\ Assume that a HF-minimizer exists with J p HF = 2V. Let 
p TF be the minimizer of the TF-energy functional of the neutral atom with nuclear 
charge Z. Then for R > to be chosen 

N= f p TF (x)dx+ / ( / u HF (x)-p TF (x))dx + f p HF (x)dx. (110) 

J\x\<R J\x\<R J\x\>R 



EXCESS CHARGE FOR PSEUDO-RELATIVISTIC ATOMS IN HARTREE-FOCK THEORY 35 



By Theorem 1 1 . 1 71 we know that there exist universal positive constants e, aa,CM 
and C$ such that for all a < ccq and x £ M. 3 

-i+e 



l*fxf(x) 



<J. TF 



(x)| <C*|x| 



C 



M ■ 



(111) 



Let Zq be such that Zqckq — k. Then a < ao corresponds to Z > Zq. Let us 
choose R such that C$i?~ 4+e = C M - Then from (fTTUf . (ITTT|) and Lemma gTJ for 
all Z > Zq we find 



N < 



p"(x) dx + 2C$i? 



-3+e 



C(l + C$i? e )(i?" 3 + 1) < Z + Q. 



|x|<_R 



The claim follows choosing Q = max{Q, + 1}. 



□ 



Proof of Theorem \1.18[ Let /9 HF be the density of the HF-minimizer in the neutral 
case TV = Z. We have 



(^(x)-^(x))dx 



|x|>fi 



(p HF (x)-p-(x))dx 



TF/ 



|x|<fl 



< 



C^R 3 



CmR, 



where in the last step we have used Theorem 11.171 Notice that for Z sufficiently 
big a < ao where «o is the constant given in Theorem ITTT71 By the TP-equation, 
Theorem 1 1.1 21 we then find 



27T 2 

3 4 — k~R 3 — C$i? 3 ^~ e — CmR^L 



|x|>_R 



p HF (x)dx<3 4 ^ir 3 - 



C^R 



-3+e 



from which the claim follows directly by the definition of HF-radius. 



-CmR, 



□ 



Proof of Theorem XTM. Since E KF (Z -l,Z)> E we (Z, Z) the ionization energy is 
bounded from below by zero. If Z is smaller than a universal constant then we can 
also bound the ionization energy with a universal constant using Theorem 12.111 

It remains to estimate from above the ionization energy when Z is larger than a 
universal constant. We first construct a density matrix 7 such that Tr[7] < Z — 1. 
Let 6- := (1 — $r(i-A)) 5 f° r r ' ^ positive parameters and 9 r defined in Definition ^. 41 
We consider the density matrix 7^ F := #_7 HF 6>_ where 7 HF is the HF-minimizcr 
in the neutral case. By an opportune choice of r we will then have TrfrP] < Z- 1. 
Indeed, 

Tr[ 7 HF ] = I P HF (x) rfx - / ^ 2 (1 _ A) (x)p HF (x) d*<Z- ( p HF (x) rfx. 

We now choose A = i. Let i? > be such that Cm = C$i? _4+e where Cm,C*,£ 
are the constants in Theorem 11.171 Then i? is a universal constant. We consider 
Z large enough so that PqZ ~a < i? where /3o is the constant in Theorem 11.121 
This gives that Z has to be larger than some universal constant. For r such that 
l\Z-* <r < i? by Theorem [LTT] we find 

|$f x f (x) - $T x F ( x )| < 2C $ |x|- 4+e for all |x| < r. 
Since J p TF = J" /0 HF , by the choice of r and Lemma I4TT1 we get 



(x) dx 



|x|>r 



> TF fxWx- 



> 



|x|>r 



|x|>t 



(/*(x)-p<"(x))dx 



|x|<r 



„TF/ 



(x) dx - 2C<s,r~ 3+e > CV~ 3 - 2C< [ ,r~ ;i+£ (112) 



3+e,- 
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In the last step we used the TF-equation, Corollary 11.131 and that r > /3qZ 3 . 
Finally, it follows from (|112[) by choosing r sufficiently small that J| x | >r P HF > 1 
and hence that Tr[ 7 ^ F ] < Z - 1. We may choose r sufficiently small by taking Z 
large enough. Notice that r can be choosen universally and so Z has to be larger 
than some universal constant. 

By the last estimate in the proof of Theorem 14.71 we find 

£ HF ( 7 HF ) < £ HF ( 7 HF ) - £ A { 1 f F ) + K, 

with TZ and 7^ F as defined in the statement of Theorem 14.71 Since £ HF (7^ F ) > 
E HF (Z- 1, Z) and £ HF ( 7 HF ) = E HF (Z, Z) it remains to prove that -£ A (^ F )+U 
is bounded from above by some universal constant. Here we use repeteadly that r is 
a universal constant. By estimate (1921 we see that 7Z < Cr~ 7 a universal constant. 
To estimate from below £ A (^ F ) we first leave out the kinetic energy term and the 
direct term since these are positive. Moreover, since $^ F is harmonic for |x| > r 
and tends to zero at infinity we see that 

$? F (x) < o sup $« F (y) < A sup <S>J F (y) + f, sup ^(y) - $« F (y)|, 

|X| |y| =J . |X| |y| =r |X| |y| =r 

which is bounded by C"/|x|, C a universal constant, by Theorem 11.171 and Corol- 
lary [TTT31 It then follows that 

£ A {lf F ) > ~ Tr[^ 7r HF ] > -- / P HF (x) dx, 
I • I r J\ xl>r 

that is bounded from below by a universal constant using Lemma 14. II □ 

Proof of Theorem \1.2(A Let ao be the constant appearing in Theorem 11.171 and Zq 
be such that aoZg = n. The claim follows directly for Z < Zq since both functions 
are bounded for |x| large, while for |x| small the functions are bounded by a constant 
times |x| _1 . 

The case Z > Zq corresponds to a < ao and for such values of a we can use the 
result in Theorem II .171 We separate the case small x, intermediate x and large x. 
Once again, comparing with the proof in the non-relativistic case ('|23j) we have to 
do an extra splitting for small x. 

By the definition of the mean field potential and Proposition 12.81 we find 

(x)-^ F (x)| < / (p TF (y) + p HF (y))( --) + 4\\p TF -P» F \\c. 

J\x-y\<s V l x -y| sJ S2 

Since p TF is bounded in La-norm, we find using Holder's inequality, Corollary 1 1.1 51 
and Lemma l3.1l that 

w 



rF (x)-^ HF (x)| < / ^fr)(-*—-l)+ C [8*zl+a-tz^*). (113) 

^|x-y|<s V |x-y| SJ 

For the integral with the HF-density we need to split the region where the HF- 
density is bounded in L 3 -norm from the one where it is bounded in is -norm. 
Proceeding as in the proof of Lemma 13.21 (from (1331) to (I3T1) replacing the integrals 
on A(|x|,fe) with integrals on |x — y| < s) using the results of Lemma T3. II we get 
with R £ (0, s) to be chosen 

I <Q HF (y)f , 1 , --) <C(zhi +R-i(aZh% +Zi +Rizi). (114) 
i|x- y |< s v |x-y| sJ 

Recall that Za — n is fixed. Choosing s such that zis^ = Z^ (i.e. s = Z~^) and 
R such that R~*Z = R^Z? (i.e R — Z~z; notice that R < s) we get from (I113[) 
and (fTTI)) 

|(^ TF (x)-^ HF (x)| < C(zi+Zi), 
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The claim follows from this inequality for x € R 3 such that |x| < {5qZ~^?- for 

7 > 0. We consider 7 < ^3. 

1+7 

If |x| > PqZ ~ then proceeding as for very small x and as in the proof of 
Theorem [3751 up to inequality (gTj]) we get for t e (^±2, f ), I > t and R < (3 Z~ l 

|<p TF (x)-^ HF (x)| < C{sr.zi + s -iz 1+ * +R-i S *Z + Z^ 3 ~ t '>). 

1 7 

Here we have also used that Za is a constant. So choosing s such that s^Zs = 
Z*(3-i) (i.e. s = Z^-l*), i? such that = z^ 3 "*) (i.e. R=Z~i+h t ) 

and optimizing in t (i.e. t = -| + we obtain 

|^ TF (x)-^ HF (x)|<CZ^i^. (115) 

Notice that t > -Mp, R < s by the choice of £ and that R satisfies the condition 
R < (3qZ~ 1 , I > t, for Z sufficiently big. The claim then follows from (|115[) for 
x e R 3 such that |x| 1+5 < p Z~i for S < We fix 5 = i^. 

We turn now to study intermediate x. Let D < 1 be such that Cm < C$D~ 4+£ 
with Cm,C§>,£ the constants in Theorem 1 1.1 71 Then for all x such that |x| < D 

|$f x F (x)-^ x F (x)|<2C 3> |xr 4 + e . 

Moreover we choose D such that Lemma \A . 1 1 1 holds . Let x be such that ftoZ~i < 
| x |i+<5 < £>t+£ w ith < ^ < S. We set r = |x| 1+ ^. Then /3 ^"3 < r < D. We 
write ^ TF (x) - <^ HF (x) = ^ TF (x) - <^ TF (x) + ^, OTF (x) - ^ HF (x) with ^ TF the 
mean field potential of the OTF-problcm defined in Subsection 14.31 By the choice 
of r and D and Lemma T4. Ill we get since |x| > r = |x| 1+M 

|^ TF (x)-^ OTF (x)|<qxr 4 -^, (116) 

for |x| > r with ( — (7 + \/73)/2. For the other two terms we see 

and proceeding as for small x with the Coulomb-norm estimate Proposition ^. 81 by 
Lemma 14.141 and inequality (|100[) 

|^ HF (x) - <^ TF (x)| < C(^r + + R-har- 7 f* + ^"V 7 + Jr^+ is 3 a 2 ) . 

Choosing v = PqT^ 1 ^ , so that va < k < 2/n, s such that s^>r~i~ = r~i + ^s~^ 
(i.e. s — r 1+ 2i), and choosing R such that the two terms where it appears are 

0_l_Q ^ — 

equal (i.e. R = r notice that R < s) we get 

|^ HF (x) - ^ r OTF (x)| < C(r- 4 +A +r- 4+3 W), 

since ar '+f is bounded and r < 1. Collecting together the inequality above and 
(TITo) and using that r = |x| 1+ ^ the claim follows for /3 Z~1? < |x| 1+<5 < . We 

fix fi = 5/2. 

1+5 

It remains to study the case of large x, i.e. |x| > D'+c with D,5,/j, universal 

1+5 

constants. For simplicity of notation we fix the universal constant A :— . We 
first notice that 

r tf 1 \ HF/ \ 

^( x) _ „TF (X) = $ HF (X) _ $ TF (X) + / g fr) ^ 

j| y |>ix| i x yi 
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The difference of the first two terms is bounded by a universal constant for |x| > A 
by the result in Theorem 11.171 To estimate the last integral we split it as follows 



|y|>|x| 



|x-y| dV ~ /|y|>|x| | x - y r y+ /|y|>|x| Ix-yl^ 

|x-y|<l |x-y|<l 



(p Lt (y) + P tit (y))dy. 

|y|>|x| 

Since |x| > A the third term on the right hand side is bounded by a universal 
constant by Lemma [4.11 (for /9 HF ) and Corollary 11.131 (for p TF ). We estimate the 
first term by Holder's inequality and Corollarv ll.151 We get a bound on the second 
term proceeding as in (|100[) (using Theorem 12. lOj) and choosing v = ^ and R = 1. 
We obtain 

/|y|>x | x -y| 

|x-y|<l 

Then there exists a universal contant A' such that |</? HF (x) — < / 9 TF (x)| < A' for 
|x|>A □ 



Appendix A. Technical lemmas 

Proof of (jTHJ) By the definition of the function G a the inequalities in (fT5]) are 
equivalent to the following ones 

ft 4 min{|t, 1} < - ft 3 < 2t 4 min{| t, 1} for t > 0. (Al) 

As before we use the substitution t — a(p/C)i . 

The estimates in (jAll) follow directly from the study of the function g separating 
the cases t < | and t>\. 

Proof of Remark 14.21 Using the estimate on K% given in ([T5|) we find 

//xeS r (ft,ft) ^(oT^x- y|) 2 dxdy 
yeS r (/3 3 ,/3 4 ) 

/"/■ e -" _1 |x-y| 

< (16) 2 a 4 // _ ^ (a q \ — n~ dxdy 

- v ' jjxeE r (ft,ft) |x-y| 4 

yeS r (/3 3 ,/?4) 

< (16) 2 a 4 e- Q_lr(/33 ^ /32 ) 47r / p- 2 dp / rfx, 

since |x — y| > (^3 — /?2)?"- The claim follows computing the two integrals. 

A.l. Fourier transform. In the present sub-section we present our notation for 
the Fourier transform (as in [ID])- Given / E L 2 (M 3 ) we denote its Fourier transform 
by 

/(P) = ^(/)(P) := -V / e ipx /(x)rfx. 

Let /, 9 e L 2 (R 3 ). The following formulas hold: 

(1) nf*9)(p) = (27r)t/(p)5(p); 

(2) 7-(/ 9 )(p) = (2 7 r)-i(/*3)(p); 

(3) if <?(x) = e- A l x l 2 then g(p) = (2A)-2 e -|pl 2 /(4A) . 

(4) |x|- Q = 7T2 (r(f))- 1 f +o ° e-^W^Af "MA for < a < n (see [H page 
130]). 
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Moreover, 



/(x)\_ i r hp) , 

¥) (k) ^lRF*' 



Appendix B. Large Z-behavior of the energy 

In [2T] the author studies the large Z-behavior of the ground state energy for 
problem JT]). In this work we are going to use the same construction in several 
points (Lemmas 13.11 l4~T2l Theorem 13.31 ....) and with, in certain cases, a slightly 
different Hamiltonian. For convenience we repeat here the main ideas of the proof. 
We do it as it is needed in the proof of Theorem 13.31 since in this case the proof 
is more involved. We remark that in our proof we use a localisation less than in 
[2T] . Thanks to Theorem [2.10l and pMl Theorem 2.8] it is sufficient to consider the 
region near the nuclei and the one far away from the nuclei. There is no need for 
an intermediate region. 

Proposition B.l. Let Za — n be fixed with < k < 2/ir and Z > 1. Let us 

consider P € R 3 , with |P| > fiZ-^ for (3 > and n € (0,4/5). Let Z > v > 
and R > be such that R < (3Z~ l /4 for some — ^ < I. Moreover, let p TF denote 
the minimizer of the TF- energy functional of a neutral atom with nucleus of charge 
Z . Consider the Hamiltonian 

N Z 1 

H P := £ (a-T( Pi ) - — - ^ (*)) + £ R _^, (B2) 

i — 1 i < j 

acting on Af =1 L 2 (R 3 ; C q ). 

Then for all t G (^^7 min{i, |}) and i/j e A*^ 2 ^ 3 ), with \\i/>\\ 2 = 1 , 

(^f p v») > ^ TF (p TF ) - co?* + /r 2 )^ 

C depending only on q and k. 

Proof. Since £ TF (p TF ) = —e$zi (see (JT3J) ) to prove the claim it is sufficient to show 
that the TF-energy gives a lower bound to the quantum energy modulo lower order 
terms. In the proof we first reduce to a one-particle operator. Then we localize 
the energy separating the contribution from the regions near the nuclei from the 
contribution from the region far away from them. Finally we study the contribution 
of each of these terms. The main contribution to the energy is given by the region 
far away from the nuclei. This region will give the TF-energy. 
In the following, s = (3 - 1)/4 (t < s < 2/3). 

In the proof C denotes a generic positive constant depending only on q and k. 

Reduction to a one-particle problem. We are going to estimate from below Hp by 
a one-particle operator. This allows us to consider only Slater determinants when 
minimizing the energy. 

Let g £ Cg°(R 3 ), g > be spherically symmetric with supp(g) C -Bi(O) and such 
that H0H2 = 1. Starting from these g we define $ s (x) := ([3/(8Z s ))- 3 g 2 (8Z s x/ f3). 
Then by Newton's theorem 



E - § J J — — dxdy 

$ a (xj -x)$ 8 (xj -y) N ff $ s (x)$ s (y) 



N 

i e 11 — -;'„T~" *■* - * IS ^f 1 *■* 
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and introducing p £ L 1 (R 3 ) n L^R^ p > 0, to be chosen 

1 r r (Eti M*» - *) - gWKgjLi <M*j - y) - gfr)) dxdy 

2 7r3 7r3 |x - y| 

^ Js? 7 R 3 |x - y| J R3 J R 3 |x - y| 

> ]Tp*$ s * — -D^J-Cll^lllJVjS- 1 ^'. (B3) 

2—1 11 

In the last inequality we use that the first term on the left hand side of f|B3|) is 
non-negative and that 



|x - y| J R 3 j K 3 |x - y| 



by definition of $ s and Hardy-Littlewood-Sobolev's inequality. Hence 

N Z 1 

- aiic? 2 !!!^- 1 ^ 8 . (B4) 



Choice of the localization. The localization will be given by the following func- 
tions xi,X2 G 



(l if|x|<^Z-*, fl if |x-P|<|/?Z-«, 

XlW := { if |x| > X2W := { if |x- P| > (B5) 

and X3 G C°°(R 3 ) such that £)? =1 X 2 (x) = 1 for all x £ R 3 . Moreover we ask that 
IIVxilU, HVxalU, llVxalU < 2 5 r 1 Z t . (B6) 

Here i is the parameter given in the statement of the proposition. Notice that 
by the assumptions on R and P the functions defined above give a well defined 
partition of unity of R 3 . Moreover, Br(P) is a subset of {x £ R 3 : X2(x) = 1}. 

The localization in the energy expectation. We insert now the localization in 
the energy expectation. As already observed, since we reduced the operator to a 
one-particle operator in the energy expectation it is sufficient to consider Slater de- 
terminants: i.e. if) = u±A- ■ -Ait at with {ui}f =1 orthonormal functions in L 2 (R 3 ,C 9 ). 
We may assume that iii £ #2(R 3 ,C«) for i = 1,...,N. 

From (|B4I) and Theorem 12. II we find with ip = u\ A • • • A itjv 

AT 3 

i=i j=i 

N 3 

-a-^^K.Lj-u,), (B7) 

i=l 3=1 

with 

»=-«- 1 r(P)-n- ! T^p+'»»-n- 

and Lj is the operator (defined in Theorem 12.11) that gives the error due to the 
localization in the kinetic energy. We first estimate this error term. Using the 
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definition of Lj we find for all j £ {1, 2, 3}, i £ {1, . . . , N} 

{uuLjUi) < ?Lgp7xj\\L J J ^ 2 (a _1 |x-y|)|u J (y)||M 4 (x)| dxdy. 
We then obtain by using Schwarz's inequality 

N 3 -3 3 N 



i=l j=l j=l i=l 

(B8) 

since from ([15} 

/ K 2 (a- 1 \z\) dz = a 3 [ K 2 (\z\) dz = Aira 3 [ t 2 K 2 (t) dt = 6ir 2 a 3 . (B9) 
Collecting together (JBT]) and flES} we get 



iV 3 

H P ^) > J^tem, fcxf«i) - - c/3- 2 z 1+2t - c/r ^ 7/4 - t/4 - (bio) 

i=l 3=1 

Here we used that N < 2Z + 1, the choice of s and that we may choose g such that 
IIVsHI < 2tt. 

Near the nuclei. When j = 1 in the summation in the first term on the right 
hand side of (|B10|) wc find 

N N „ 

i=l »=1 ' ' ' 

since Xb h (P)Xi = by the choice of Xi) an d the term $ s * p * A- is non- negative. 
Then by Theorem 1 2 . 1 01 we find 

N Z 
/^(Xi^hxiUi) > Tr[a _1 r(p) - rTX|x|<i0Z-*]- 

i=l ' ' ' 

> -C*/3 1/2 Z 5/2 - t/2 -CK 2 Z 2 . (Bll) 

To estimate from below the term corresponding to j = 2 in the sum on the 
right hand side of (|B10[) we use PU Theorem 2.8]. Here we need the result in (53] 
(instead of Theorem I2.10j) because of the presence of the two nuclei. Notice that 
Theorem 12.101 can be extended to include also different nuclei. We have 

N N 
^2(X2Ui, h X 2Ui) > ^(X2"i, (a _1 T(p) - r— - —— — XB R (P))X2"i) 
i=l i=l ' X ' ' X I 

> Trla-^Cp) - j^|X|x-p|<^z-« - [x z p] Xa fi (P)]-, 
and by [MJ Theorem 2.8] we get 

N r / 75/2 ^4 \ 

^ 7i/3z-'>|x-p|>a Vl x l 5/2 l x rv 

/" / 3 ^ \ 

- c I>,x-P,>AR^F + a 'FpFj dx 

> -CV 1/2 Z 2 - C(3 1/2 Z 5/2 ~ t/2 _ Ck 2 Z 2 . (B12) 
Here we used that t < I and Za = k. 
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The outer zone. This region gives the main contribution to the energy. The term 
in (|B10j) that we still have to study is 

N 

E^ 3 " 1 ' h X 3Ui) - D(p) (B13) 
i=i 

We start by estimating the first term in (|B13|) using coherent states. 

We consider again the function g g C£°(R 3 ) introduced at the beginning of the 
proof and we define the function 

<? s (x) := ((3/(8Z s )yig(8Z s x/(3) = <&f (x), (B14) 

with s the same parameter as before. For simplicity of notation we write V :— 
Z/\x\ - p* l/|x|. Then 

Z 1 1 Z 

— - p * $ s * — =v*$ s - Z<$> s * — + — . 
|x| |x| |x| |x| 

Since supp(<7 s ) fl supp(x3) = by Newton's Theorem we find 

N N 

X)(*3«<, hxiUi) = X)(X3«i, (a _1 T(p) - V * $ s ) X 3Ui). (B15) 
i=l i=l 

We consider the coherent states gf ,cl defined for p, q G R 3 by 

3 ^(x)= gs (x-q) e -' ; P x . 
The following formulas hold for / e H^(R 3 ,C) 



(/,/) 



1 / dp / dq(/ >5 ?*) (</?*/), 

V ; JR 3 JR 3 



and 



(/,^*ff s 2 /) = J2W f dp f dqV(q)(/,3?' q )(sTV) (B16) 

JR 3 JR 3 

(/,T( P )/) = i / dp/ dqr( P ) 

JR 3 JR 3 

/ dx / dq7^(£,/)(x), (B17) 

JR 3 JR 3 



where L g has integral kernel 



-2 



. 2 K 2 {a- l \*-y\) 



- 4?r2 HJ ysw HJ\ | x -y| 2 

Using these formulas we can rewrite (1B15[) as follows 

N 
i=l 

q N 

EE 

3 = 1 i=l 



= f dp f d q (T(p)-aV( q ))^2Y,\(x 3 ul,gf^\ 2 

-cT 1 V / dx / dqx3U 4 (x)(£ q X3" i )(x), (B18) 

i=1 JR 3 JR 3 



Here u\ is the j-th spin component of Ui. We start by estimating the error term, 
the last term on the right hand side of (|B18[) . From the definition of L q it follows 

£«(x,y) < ^2-||V5 s ||^i^ 2 (a _1 |x-y|)(Xsupp( ffs )(x- q) + Xsu PP ( 9s ) (y - q)), 
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and by the definition of the function g s 

f L 9 (x,y)dq<q|V.g||^a- 2 r 2 ^ 2 ^ 2 (a- 1 |x-y|). 

JV. 3 

By the estimate above, Schwarz's inequality, (|B9I) and the choice of s we find 



a ^}Z / dx / dcixzu^iL^u^KCWVgWl^Z^-^N. (B19) 



It remains to study the first term on the right hand side of (|B18I) . In order to 
get an estimate from below we consider only the negative part of the integrand. 
Moreover, since if |q| < f3Z^ t /8 then supp(x3(? P ' q ) = (because Z~ l > Z~ s since 
s > t) we find 



r r q N 

j^a- 1 dp dq(T(p)- a y(q))££|( X3 «' q ) 



> 7^ a dc ± dp(T(p)-aV(q)) = ...,(B20) 



where we also use that $D j=1 l(X3 u ii5!r q )| 2 < 1 (Bessel's inequality). We split now 
the integral as a sum of two terms 

- = * a " 1 //iW-^<o* rt, ™- a ^ 

|q|>^z-' 

+ (2^ / / f |p| 2 >^( q) >T(p) ^ ( T (P) - «^). (B21) 

|q|>A/92-* 

We consider these two terms separately. The second term in (|B21|) gives a lower 
order contribution. Indeed 

(sfF^ 1 / / f |p|^> Q y (q )>T(p) d * d P ( T (P) ~ 
|q|>s/3z-' 



^ (2-) 3 yy(c t 2 [y( q )] 2 ++ 2[V-(q)] + )i>|p|>(2[y(q)] + )i ^ - • • • ' 

|q|>|/3z-* 

and computing the p-intcgral 

■ =-C [ rfq [F(q)]|((l + ^[?(q)]+)* - 1) = . . . . 

Using (1 + ar)i < 1 + |x + |x 2 and that [V"(q)] + < Z/|q| we get computing the 
integral 

••• = -Ca a Jl q |> ii9 z- t dq[V'(q3]i(l + ^[^(q)] + ) 

(BzzJ 

> -C/3-5 K 2Z 3 / 2 +*/2 _ Z l/2+3t/2_ 

Here we use that Za = k. 



Since y/T+x > 1 + x/2 — x /8 for all cc > 0, we have 

T(p)>ai|p| 2 -a 3 i|p| 4 , 
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and, for the first term on the right hand side of (|B21I) . we obtain 

|q|>|/3Z"* 



(2jr)3 yy3ipi 2 -^(q)<° V21 1 8 

|q|>J^-* 

Computing now the integral with respect to p, we find 



— 2 2<j 
_ 157T 2 



/ [F(q)]|dq-CV* 2 / [F(q)]f<iq. (B23) 

,/|q|>i/3Z-* J|q|>i^Z-* 



We see that the second term on the right hand side of (|B23I) gives a lower order 
contribution since it is of the same order as the one in (|B22p . 

Collecting together ([Biol , (|BTT|) . (1B12)) . (|Bl~5l) . (|Bl8)) . ([Hl9]> . (|B22)) and (|B23)) 

(^ffp^> > -C(/3U/?- 2 )^ 5/2 - t/2 -ill / [^(q)] I dq-£>(p).(B24) 

JR 3 

Here we used also that A*" < 2Z + 1, the choice of s and that i < 3/5. 

Now we choose p = p TF the minimizer of the TF-energy functional of a neutral 
atom with Coulomb potential and nuclear charge Z. Hence p TF satisfies the TF- 
equation 

i(^)V F (x)t = [y(x)] +I 

since V is the TF-mean field potential. Notice that here we use that the chemical 
potential of a neutral atom is zero. By the choice of p from the TF-equation it 
follows from (jB24f that 

{<l>,H P i/,) > -C^+r^-^ + U^) 2 " I rfxp TF (x)§ 

JR 3 

-z [ ^!W dx + ^ TF ) 

= £ TF (p TF )-C(^ + f3- 2 )Z 5/2 ~ t/2 . 
The claim follows. □ 

Proposition B.2. Let p TF be the minimizer of the TF-energy functional of a 
neutral atom with nuclear charge Z. Let Za = k be fixed with < n < 2 /it and 
Z > 1. 

Then there is a constant depending only on k and q such that for all {ui}fL 1 C 
H^(R 3 ;C q ) orthonormal in L 2 (R 3 ) we have 

N 

5>„ (a-'Tfp) - ^ F )^) - D(p TF ) > £ TF ( P TF ) - CZ*+i , 

i=l 

with D(-) = D(-, •) the Coulomb scalar product. 

Proof. Since £ TF (p TF ) = —e^Zi (see lfT2"|) ) to prove the claim it is sufficient to 
show that the TF-energy gives a lower bound to the quantum energy modulo lower 
order terms. In the proof we localize the energy separating the contribution from 
the region near the nucleus to the one far away. The region far away from the nuclei 
will give the TF-energy. 

In the proof C denotes a generic universal positive constant. 
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Choice of the localization. The localization will be given by the functions Xi S 
(^(R 3 ) and X2 € C°°(M 3 ) such that: < xi,X2 < 1, Xi + xl = 1 in K 3 , 

if |x| > 3Z- 3 / 5 . (B25) 

Moreover we ask that 

||Vxi||oo,||Vxa||oo <2 2 Z 3 / 5 . (B26) 

The localization in the energy expectation. We insert now the localization in the 
energy expectation. From Theorem 12. II we find 

N 

(a' 1 r(p) - ^ F ) Ui ) - D(p TF ) (B27) 

i=l 

N 2 AT 2 

> E («" lT (p) - ^ TF )x^i) - d (p tf ) - « _1 EEk^), 

2—1 J — 1 Z— 1 j — 1 

with is the operator (defined in Theorem I2.1[) that gives the error due to the 
localization in the kinetic energy. We first estimate this error term. Since N < 
2Z + 1 we find as in (TbSI) that 



iV 2 

a' 1 J2(ui,L jUi ) < CZ e/5 N < CZ 2+1/5 . (B28) 

»=1 3=1 



iVear £/ie nucleus. Since 

JV 



i=l 

by Theorem HTU] with R = 3Z~ 3 / 5 we find 

N 

J2(xm, (a- 1 T(p) - ^ F )xiui) > -CZ 2+1 ' 5 - Cn 2 Z 2 . (B29) 



Here we use that Za = k. 

The outer zone. This region gives the main contribution to the energy. 

Let g € C5°(M 3 ), g > be spherically symmetric with supp(g) C £?i(0) and such 
that ||<7|| 2 = 1. Starting from these g we define $z(x) := (Z~ 3//5 )~ 3 g 2 (xZ 3 / 5 ) and 

ffz (x) := (Z- 3 / 5 )"i 5 (xZ 3 / 5 ) = $|(x). 
Since supp(gz) H supp(x2) = by Newton's Theorem we find 

N N 

Y J (X2U t , (a- x T(p) - i P TF ) X 2U i ) = X)(X2«i, (a _1 T(p) - f T¥ * ®z)X2Ui). (B30) 

i=i 1=1 

We consider the coherent states g^' q defined for p, q £ R 3 by 

^ q (x)=. 9z (x-q)e-^ x . 
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Using formulas (|B16I) and (IB17[) we can rewrite (|B30[) as follows 

N 

^( X 2Ui, (a-'Tip) - ^ TF * g 2 z ) X 2U t ) 
i=i 

r r q N 

= dp dq(T(p)-a^ TF (q))^^|(x 2 ^,ffni 2 

Jr 3 J r3 ]=1i=i 

-a" 1 V / dx dqx 2 u i (x)(LqX2Mi)(x), (B31) 
i=1 Jr 3 Jr 3 

Here u\ is the j-th spin component of Ui. We start by estimating the error term, 



the last term on the right hand side of (|B31|) . We find as in (|B19j) that 

/ dx I rfqAWx)(i q X2^)(x)<C||V. 9 ||L^ 6/5 ^. (B32) 

■ =1 JR 3 JR 3 

It remains to study the first term on the right hand side of (|B31|) . In order to 
get an estimate from below we consider only the negative part of the integrand. 
Moreover, since if |q| < Z~ 3 / 5 then supp(x2<?z q ) = we find 



r r q N 

rtT 1 / dp dq(T(p)-^ TF (q))^^|( X2 ui,^ q )| 

Jr 3 Jr 3 



3 = 1 i=l 



> T2^ a ^ dq dp (T(p) - a^ TF (q)) = . . .(B33) 

J|q|>Z- 3 / 5 JT{p)-a V TF (q)<0 

where we also use that J2f =1 |(X3 u i7 5z q )l 2 — 1 (Bessel's inequality). We split now 
the integral as a sum of two terms 

ill frr/\ „,TF 



i^7 a I k\ P ^-^M<o d( idp(T(p)-^ (q)) 

\q\>Z~ 3/5 

«_^v-l // rlrtrln (Tlr\\ — min^^ I 



+ ^ yy f i P iw p (c 1 )>T(p)^p( T (p)-^ ( B34 ) 

We consider these two terms separately. The second term in (|B34[) gives a lower 
order contribution. Indeed 



W7 a 1 I l^> a ^(c l )>T( P )d t idp(T(p)-a l p TF ( q )) 

|q|>Z~ 3/5 



- -(^Fy J /(^[,-] 2 ++2 [,-] + )i>iPi>( 2 [,-( q )] + )i dqdp [ ^ TF(q)]+ 

|q|>Z~ 3/5 

and computing the integral in p 



' • • = ~C [ dq ^ TF (q)]|((l + |V F (q)] + )i " 1) 

J\q\>Z~ 3 / 5 Z 

1 

integral 



'|q|>-Z- 3 / 5 

Using (1 + x) 'i < 1 + |a; + |a; 2 and that [<£ TF (q)] + < Z/\q\ we get computing the 



= -Co?\ dq[^ TF ]|(l + ^V F (q)] + ) 

J \q\>z-*/s » (B35) 



> -C K 2 Z 2 -i - Cn 4 Zi. 
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Since ^/l + x > 1 + x/2 — x 3 /8 for all x > 0, we have 

T(p)>a±|p| 2 -a 3 ||p| 4 , 
and, for the first term on the right hand side of (|B34|) . we obtain 

J hrf-^(,)<o dqdp (T(P) " atpTF{q)) * 

| q |> z -3/5 

^ i^f /4p|^-( q )<o^p(ilpl 2 -Xlpl 4 -^ TP W) 



|q|>Z 



-3/5 



Computing now the integral with respect to p, we find 

■■• = / ^ TF (q)]}rfq-Ca 2 / [^ TF (q)]J dq. (B36) 

J|q|>Z-3/5 J|q|>Z-3/ 5 

We see that the second term on the right hand side of (|B36[) gives a lower order 



contribution since it is of the same order as the one in (|B35|) . 

Starting from (|B27j) . by (|B28|) . (tB29|) . (fB32l) . f[B35|) and (|B36j) wc find 

N 

(a-^Cp) - ^ TF K) - £(p TF ) (B37) 

i=l 

> -C(Z 2+1 /5 + Z 2 + Z 2 - 1 '^ + Z 7 / 5 ) - ^ / b TF (q)]|dq-i?(p TF ). 
The result follows from the TF-equation. □ 
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